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Abstract 

This chapter aims to present an introduction to current research on the 
nature of the cosmological dark matter and the origin of galaxies and large 
scale structure within the standard theoretical framework: gravitational 
collapse of fluctuations as the origin of structure in the expanding universe. 
General relativistic cosmology is summarized, and the data on the basic 
cosmological parameters {to and Hq = 100/ikms~^Mpc~^, JIq; ^^A and 
rib) a-re reviewed. Various particle physics candidates for hot, warm, and 
cold dark matter are briefly reviewed, together with current constraints 
and experiments that could detect or eliminate them. Also included is 
a very brief summary of the theory of cosmic defects, and a somewhat 
more extended exposition of the idea of cosmological inflation with a 
summary of some current models of inflation. The remainder is a discussion 
of observational constraints on cosmological model building, emphasizing 
models in which most of the dark matter is cold and the primordial 
fluctuations are the sort predicted by inflation. It is argued that the 
simplest models that have a hope of working are Cold Dark Matter with 
a cosmological constant (ACDM) if the Hubble parameter is high (h ^ 0.7), 
and Cold + Hot Dark Matter (CHDM) if the Hubble parameter and age 
permit an = 1 cosmology, as seems plausible in light of the data from 
the Hipparcos astrometric satellite. The most attractive variants of these 
models and the critical tests for each are discussed. 



To be published as Chapter 1 of Formation of Structure in the Universe, 
Proceedings of the Jerusalem Winter School 1996, edited by A. Dekel and 
J. P. Ostriker (Cambridge University Press). 
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1.1 Introduction 

The standard theory of cosmology is the Hot Big Bang, according to which 
the early universe was hot, dense, very nearly homogeneous, and expanding 
adiabatically according to the laws of general relativity (GR). This theory 
nicely accounts for the cosmic background radiation, and is at least roughly 
consistent with the abundances of the lightest nuclides. It is probably even 
true, as far as it goes; at least, I will assume so here. But as a fundamental 
theory of cosmology, the standard theory is seriously incomplete. One way 
of putting this is to say that it describes the middle of the story, but leaves 
us guessing about both the beginning and the end. 

Galaxies and clusters of galaxies are the largest bound systems, and the 
filamentary or wall-like superclusters and the voids between them are the 
largest scale structures visible in the universe, but their origins are not yet 
entirely understood. Moreover, within the framework of the standard theory 
of gravity, there is compelling observational evidence that most of the mass 
detected gravitationally in galaxies and clusters, and especially on larger 
scales, is "dark" — that is, visible neither in absorption nor emission of any 
frequency of electromagnetic radiation. But we still do not know what this 
dark matter is. 

Explaining the rich variety and correlations of galaxy and cluster 
morphology will require filling in much more of the history of the universe: 

• Beginnings, in order to understand the origin of the fluctuations that 
eventually collapse gravitationally to form galaxies and large scale 
structure. This is a mystery in the standard hot big bang universe, because 
the matter that comprises a typical galaxy, for example, first came into 
causal contact about a year after the Big Bang. It is hard to see how 
galaxy-size fluctuations could have formed after that, but even harder to 
see how they could have formed earlier. The best solution to this problem 
yet discovered, and the one emphasized here, is cosmic inflation. The main 
alternative, discussed in less detail here, is cosmic topological defects. 

• Denouement, since even given appropriate initial fluctuations, we are far 
from understanding the evolution of galaxies, clusters, and large scale 
structure — or even the origins of stars and the stellar initial mass 
function. 

• And the dark matter is probably the key to unravelling the plot since it 
appears to be gravitationally dominant on all scales larger than the cores 
of galaxies. The dark matter is therefore crucial for understanding the 
evolution and present structure of galaxies, clusters, superclusters and 
voids. 
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The present chapter (updating Primack 1987-88, 1993, 1995-97) 
concentrates on the period after the first three minutes, during which the 
universe expands by a factor of ~ 10*^ to its present size, and ah the observed 
structures form. This is now an area undergoing intense development in 
astrophysics, both observationahy and theoreticahy. It is hkely that the 
present decade wih see the construction at last of a fundamental theory of 
cosmology, with perhaps profound implications for particle physics — and 
possibly even for broader areas of modern culture. 

The current controversy over the amount of matter in the universe 
will be emphasized, discussing especially the two leading alternatives: a 
critical-density universe, i.e. with = po/pc = 1 (see Table 1.1), vs. a 
low-density universe having 0,q ~ 0.3 with a positive cosmological constant 
A > such that Oa ^ A/{3Hq) = 1 — supplying the additional density 
required for the flatness predicted by the simplest inflationary models. (The 
significance of the cosmological parameters Oq, Hq, to, and A is discussed in 
§ ^^■) $7 = 1 requires that the expansion rate of the universe, the Hubble 
parameter Hq = 100/ikms~^Mpc~^ = 5O/150 kms~^Mpc~^, be relatively 
low, h ^ 0.6, in order that the age of the universe tg be as large as the 
minimum estimate of the age of the stars in the oldest globular clusters. If 
the expansion rate turns out to be larger than this, we will see that GR then 
requires that f^o < 1, with a positive cosmological constant giving a larger 
age for any value of f^o- 

Although this chapter will concentrate on the implications of CDM and 
alternative theories of dark matter for the development of galaxies and 
large scale structure in the relatively "recent" universe, we can hardly 
avoid recalling some of the earlier parts of the story. Inflation or cosmic 
defects will be important in this context for the nearly constant curvature 
(near-"Zerdovich") spectrum of primordial fluctuations and as plausible 
solutions to the problem of generating these large scale fluctuations without 
violating causality; and primordial nucleosynthesis will be important as a 
source of information on the amount of ordinary ( "baryonic" ) matter in the 
universe. The fact that the observational lower bound on — namely 
0.3 ^ — exceeds the most conservative upper limit on baryonic mass 
~ 0.03/i~^ from Big Bang Nucleosynthesis (Copi, Schramm, & Turner 
1995; cf. Hata et al. 1995) is the main evidence that there must be such 
nonbaryonic dark matter particles. 

Of special concern will be evidence and arguments bearing on the 
astrophysical properties of the dark matter, which can also help to constrain 
possible particle physics candidates. The most popular of these are few-eV 
neutrinos (the "hot" dark matter candidate), heavy stable particles such as 
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~ 100 GeV photinos (or whatever neutralino is the hghtest supersymmetric 
partner particle) or 10~® — 10~^ cV "invisible" axions (these remain the 
favorite "cold" dark matter candidates), and various more exotic ideas such 
as keV gravitinos ("warm" dark matter) or primordial black holes (BH). 

Here we are using the usual astrophysical classification of the dark matter 
candidates into hot, warm, or cold, depending on their thermal velocity 
in the early universe. Hot dark matter, such as few-eV neutrinos, is still 
relativistic when galaxy-size masses (~ IO^^Mq) are first encompassed 
within the horizon. Warm dark matter is just becoming nonrelativistic then. 
Cold dark matter, such as axions or massive photinos, is nonrelativistic 
when even globular cluster masses (~ 10^ Mq) come within the horizon. 
As a consequence, fluctuations on galaxy scales are wiped out by the "free 
streaming" of the hot dark matter particles which are moving at nearly the 
speed of light. But galaxy-size fluctuations are preserved with warm dark 
matter, and all cosmologically relevant fluctuations survive in a universe 
dominated by the sluggishly moving cold dark matter. 

The first possibility for nonbaryonic dark matter that was examined in 
detail was massive neutrinos, assumed to have mass 25 eV — both 
because that mass corresponds to closure density for h ^ 0.5, and because 
in the late 1970s the Moscow tritium /3-decay experiment provided evidence 
(subsequently contradicted by other experiments) that the electron neutrino 
has that mass. Although this picture leads to superclusters and voids 
of roughly the size seen, superclusters are the first structures to collapse 
in this theory since smaller size fluctuations do not survive. The theory 
foundered on this point, however, since galaxies are almost certainly older 
than superclusters. The standard (adiabatic) form of this theory has recently 
been ruled out by the COBE data: if the amplitude of the fluctuation 
spectrum is small enough for consistency with the COBE fluctuations, 
superclusters would just be beginning to form at the present epoch, and 
hardly any smaller-scale structures, including galaxies, could have formed 
by the present epoch. 

A currently popular possibility is that the dark matter is cold. After 
Peebles (1982), we were among those who first proposed and worked out 
the consequences of the Cold Dark Matter (CDM) model (Primack & 
Blumenthal 1983, 1984; Blumenthal et al. 1984). Its virtues include an 
account of galaxy and cluster formation that at first sight appeared to 
be very attractive. Its defects took longer to uncover, partly because 
uncertainty about how to normalize the CDM fluctuation amplitude allowed 
for a certain amount of fudging, at least until COBE measured the 
fluctuation amplitude. The most serious problem with CDM is probably 
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Table 1.1. Physical Constants for Cosmology 



speed of light 
solar mass 
solar luminosity 

Planck's const. 
Planck mass 



parsec 

Newton's const. 
Hubble parameter 
Hubble time 
Hubble radius 
critical density 




rape 



c 



Pc 



h 



© 



3.09 X 10^8 cm = 3.26 light years (lyr) 

6.67 X 10"^ dyne cm^ g-^ 

lOO/i km s-i Mpc-i , l/2</i<l 

h-^ 9.78 Gyr 

cH-^ =3.00/1-1 Qpp 

3H^/8TrG = 1.88 x IQ-^^h^ g cm'^ 

10.5 keV cm'^ = 2.78 x lO^/i^ Mq Mpc'^ 

3.00 X 10^° cm s-i = 306 Mpc Gyr"! 

1.99 X 10^3 g 

3.85 X 10^3 erg s'^ 

1.05 X 10-27 erg s = 6.58 x IQ-^^ s 
ihclGfl'^ = 2.18 X 10-5 g = 1.22 X 10^^ GeV 



the mismatch between supcrclustcr-scale and galaxy-scale structures and 
velocities, which suggests that the CDM fluctuation spectrum is not quite 
the right shape — which can perhaps be remedied if the dark matter content 
is a mixture of hot and cold, or if there is less than a critical density of cold 
dark matter. 

The basic theoretical framework for cosmology is reviewed first, followed 
by a discussion of the current knowledge about the fundamental cosmological 
parameters. 

Table 1.1 lists the values of the most important physical constants used 

in this chapter (cf. Barnctt et al. 1996). The distance to distant galaxies 
is deduced from their rcdshifts; consequently, the parameter h appears in 
many formulas where the distance matters. 



It is assumed here that Einstein's general theory of relativity (GR) 
accurately describes gravity. Although it is important to appreciate that 
there is no observational confirmation of this on scales larger than about 1 
Mpc, the tests of GR on smaller scales are becoming increasingly precise, 
especially with pulsars in binary star systems (Will 1981, 1986, 1990; Taylor 
1994). On galaxy and cluster scales, the general agreement between the mass 
estimated by velocity measurements and by gravitational lensing provides 
evidence supporting standard gravity. There are two other reasons most 
cosmologists believe in GR: it is conceptually so beautifully simple that it is 
hard to believe it could be wrong, and anyway it has no serious theoretical 
competition. Nevertheless, since a straightforward interpretation of the 
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available data in the context of the standard theory of gravity leads to the 
disquieting conclusion that most of the matter in the universe is dark, there 
have been suggestions that perhaps our theory of gravity is inadequate on 
large scales. They are mentioned briefly at the end of this section. 

The "Copernican" or "cosmological" principle is logically independent of 
our theory of gravity, so it is appropriate to state it before discussing GR 
further. First, some definitions are necessary: 

• A co-moving observer is at rest and unaccelerated with respect to nearby 
material (in practice, with respect to the center of mass of galaxies within, 
say, 100 h-^ Mpc). 

• The universe is homogeneous if all co-moving observers see identical 
properties. 

• The universe is isotropic if all co-moving observers see no preferred 
direction. 

The cosmological principle asserts that the universe is homogeneous and 
isotropic on large scales. (Isotropy about at least three points actually 
implies homogeneity, but the counterexample of a cylinder shows that the 
reverse is not true.) In reality, the matter distribution in the universe is 
exceedingly inhomogeneous on small scales, and increasingly homogeneous 
on scales approaching the entire horizon. The cosmological principle is in 
practice the assumption that for cosmological purposes we can neglect this 
inhomogeneity, or treat it perturbatively. This has now been put on an 
improved basis, based on the observed isotropy of the cosmic background 
radiation and the (partially testable) Copernican assumption that other 
observers also see a nearly homogeneous CBR. The "COBE-Copernicus" 
theorem (Stoeger, Maartens, Sz Ellis 1995; Maartens, Ellis, & Stoeger 1995; 
reviewed by Ellis 1996) asserts that if all comoving observers measure the 
cosmic microwave background radiation to be almost isotropic in a region 
of the expanding universe, then the universe is locally almost spatially 
homogeneous and isotropic in that region. 

The great advantage of assuming homogeneity is that our own cosmic 
neighborhood becomes representative of the whole universe, and the range 
of cosmological models to be considered is also enormously reduced. The 
cosmological principle also implies the existence of a universal cosmic time, 
since all observers see the same sequence of cosmic events with which to 
synchronize their clocks. (This assumption is sometimes explicitly included 
in the statement of the cosmological principle; e.g., Rindler (1977), p. 203.) 
In particular, they can all start their clocks with the Big Bang. 

Astronomers observe that the redshift z = (A— Ao)/Ao of distant galaxies is 
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proportional to their distance. We assume, for lack of any viable alternative 
explanation, that this redshift is due to the expansion of the universe. 
Recent evidence for this includes higher CBR temperature at higher redshift 
(Songaila et al. 1994b) and time dilation of high-redshift Type la supernovae 
(Goldhaber et al. 1996). The cosmological principle then implies (see, for 
example, Rowan-Robinson 1981, §4.3) that the expansion is homogeneous: 
r = a{t)ro, which immediately implies Hubble's law: v = r = aa~^r = Hqt. 
Here tq is the present distance of some distant galaxy (the subscript "0" in 
cosmology denotes the present era), r is its distance as a function of time 
and V is its velocity, and a{t) is the scale factor of the expansion (scaled 
to be unity at the present: a(to) = !)• The scale factor is related to the 
redshift by a = (l + z)"^. Hubble's "constant" H{t) (constant in space, but 
a function of time except in an empty universe) is H{t) = da~^. 

Finally, it can be shown (see, e.g., Weinberg 1972, Rindler 1977) 
that the most general metric satisfying the cosmological principle is the 
Robertson- Walker metric 



where the curvature constant k, by a suitable choice of units for r, has the 
value 1,0, or -1, depending on whether the universe is closed, flat, or open, 
respectively. For k = 1 the spatial universe can be regarded as the surface 
of a sphere of radius a{t) in four-dimensional Euclidean space; and although 
for A; = or —1 no such simple geometric interpretation is possible, a{t) still 
sets the scale of the geometry of space. 

Formally, GR consists of the assumption of the Equivalence Principle (or 
the Principle of General Covariance) together with Einstein's field equations, 
labeled (E) in Table 1.2, where the key equations have been collected. The 
Equivalence Principle implies that spacetime is locally Minkowskian and 
globally (pseudo-)Riemannian, and the field equations specify precisely how 
spacetime responds to its contents. The essential physical idea underlying 
GR is that spacetime is not just an arena, but rather an active participant 
in the dynamics, as summarized by John Wheeler: "Matter tells space how 
to curve, curved space tells matter how to move." 

Comoving coordinates are coordinates with respect to which comoving 
observers are at rest. A comoving coordinate system expands with 
the Hubble expansion. It is convenient to specify linear dimensions in 
comoving coordinates scaled to the present; for example, if we say that 
two objects were 1 Mpc apart in comoving coordinates at a redshift of 
z = 9, their actual distance then was 0.1 Mpc. In a non-empty universe 



ds^ = c^dt^ - a{tf 



dr"^ 



+ r^ {sin^ 9d(j)^ + de^) , 



(1.1) 



1 — kr"^ 
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GR: Matter tells space 

HOW TO CURVE, 



Curved space tells 
matter how to move. 



(E) Ri^" - ^Rg>"' = -8ttG Tf" - Agi^" 



COBE - Copernicus Th: If all observers measure nearly isotropic CBR, then 
universe is locally nearly homogeneous and isotropic - i.e., nearly FRW. 



PRWE(0«, 5 = |0.44 

FRW E(m) — + r_ = -SttGp - ^ + a 
a a'^ a^ 



Ho = 100/ikms"^Mpc-2 
= 50/i5okms~^Mpc~^ 

^ l = no-^+flAwithHo = ^, ao = l, no=f^,^A=^, 



Hi 



SttG 



: 0.70 X 10"/iioMQMpc" 



E{ii) - E{m) 



— = -— Gp-87rGp+-A 
a 3 3 



1^0 



Divide by 2^;(00) ^ = - [- = ^-^ - 



a a' 



' a 



Stt^ k A 
—Gp - ^ + - 



a3 



tQ = Ho^f {Qo,^a) = 9.78/i-iGyr 



k 


f f^A 


/(i,o) 


2 

~ 3 


/(0,0) 


= 1 


/(o,i) 


= OO 



d 

[E{00)a^]' vs. E{ii) =^ g^^P"'^'^ = ("continuity") 

Given eq. of state p = p{p), integrate to determine p{a), 

integrate i?(00) to determine a{t) 

Examples: P ~ =^ p = poa^^ (assumed above in go, h eqs.) 
p = ^, k = ^ p (X a~'^ 



with vanishing cosmological constant, the case first studied in detail by 
the Russian cosmologist Alexander Friedmann in 1922-24, gravitational 
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attraction ensures that the expansion rate is always decreasing. As a result, 
the Hubble radius Ruit) = cH{t)~^ is increasing. The Hubble radius of a 
non-empty Friedmann universe expands even in comoving coordinates. Our 
backward lightcone encompasses more of the universe as time goes on. 



1.2.1 Friedmann- Roberts on- Walker Universes 

For a homogeneous and isotropic fluid of density p and pressure p in 
a homogeneous universe with curvature k and cosmological constant A, 
Einstein's system of partial differential equations reduces to the two ordinary 
differential equations labeled in Table 1.2 FRW E(00) and E(ii), for the 
diagonal time and spatial components (see, e.g., Rindler 1977, §9.9). 
Dividing E(00) by Hq, and subtracting E(00) from E(ii) puts these equations 
into more familiar forms. Dividing the latter by 2E(00) and evaluating all 
expressions at the present epoch then gives the familiar expression for the 
deceleration parameter go in terms of 0,q and 

Multiplying E(00) by a^, differentiating with respect to a, and comparing 
with E(ii) gives the equation of continuity. Given an equation of state p = 
p(p), this equation can be integrated to determine p{a); then E(00) can be 
integrated to determine a{t). 

Consider, for example, the case of vanishing pressure p = 0, which is 
presumably an excellent approximation for the present universe since the 
contribution of radiation and massless neutrinos (both having p = pc^/S) 
to the mass-energy density is at the present epoch much less than that of 
nonrelativistic matter (for which p is negligible). The continuity equation 
reduces to (47r/3)pa'^ = M = constant, and E(00) yields Friedmann's 
equation 

= kc^ + . 1.2 

a 6 

This gives an expression for the age of the universe to which can be integrated 
in general in terms of elliptic functions, and for A = or /c = in terms of 
elementary functions (cf. standard textbooks, e.g. Peebles 1993, §13, and 
Felton & Isaacman 1986). 

Figure (a) plots the evolution of the scale factor a for three interesting 
examples: (J1o,^^a) = (1,0), (0.3,0), and (0.3,0.7). Figure |L1| (b) shows how 
t/tn depends on both for A = (dashed) and f^A = 1 — f^o (solid). 
Notice that for A = 0, to/tn is somewhat greater for = 0.3 (0.81) 
than for = 1 (2/3), while for Qo = 1 ~ = 0.3 it is substantially 
greater: to/^H = 0.96. In the latter case, the competition between the 
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Fig. 1.1. (a) Evolution of the scale factor a{t) plotted vs. the time after the present 
(t — ta) in units of Hubble time tn = H^^ = 9.78/i~^ Gyr for three different 
cosmologies: Einstein-de Sitter (Qq — l,f2A — dotted curve), negative curvature 
(f^o = O.SjfiA = 0: dashed curve), and Iow-JIq flat {flo = 0.3, fl^ = 0.7: solid 
curve), (b) Age of the universe today to in units of Hubble time tn as a function 
of fio for A = (dashed curve) and flat + ^Ia — 1 (solid curve) cosmologies. 

attraction of the matter and the repulsion of space by space represented 
by the cosmological constant results in a slowing of the expansion at 
a ~ 0.5; the cosmological constant subsequently dominates, resulting in 
an accelerated expansion (negative deceleration qo = —0.55 at the present 
epoch), corresponding to an inflationary universe. In addition to increasing 
to, this behavior has observational implications that we will explore in 
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1.2.2 Is the Gravitational Force oc r at Large r? 

Back to the question whether our conventional theory of gravity is 
trustworthy on large scales. The reason for raising this question is that 
interpreting modern observations within the context of the standard theory 
leads to the conclusion that at least 90% of the matter in the universe is dark. 
Moreover, there is no observational confirmation that the gravitational force 
falls as on large scales. 

Tohline (1983) pointed out that a modified gravitational force law, with 
the gravitational acceleration given by a' = (GM/u^/r^) (1 + r/d) , could 
be an alternative to dark matter galactic halos as an explanation of the 
constant- velocity rotation curves of spiral galaxies. The mass is written as 
^lurn to emphasize that there is not supposed to be any dark matter.) 
Indeed, this equation implies v"^ = GAIin„i/d = constant for r ^ d. 
However, with the distance scale d where the force shifts from to 
taken to be a physical constant, the same for all galaxies, this implies 
that Miujn oc v'^, whereas observationally M;„^ oa Lb oc v"' with a ~ 4 
( "Tully-Fisher law"), where Lb is the galaxy luminosity in the blue band. 

Milgrom (1983, 1994, 1995; cf. Mannheim & Kazanas 1994) proposed 
an alternative idea, that the separation between the classical and modified 
regimes is determined by the value of the gravitational acceleration a' rather 
than the distance scale r. Specifically, Milgrom proposed that 

a' = GMiumr'^, a' > a'^, a'^ = GMiumr'^a'o, a < a[) (1.3) 

where the value of the critical acceleration Oq ~ 8 x 10^^/i^cms~^ is 
determined for large spiral galaxies with Mi^^ ~ 10^^ M©. (This value for 
Oq happens to be numerically approximately equal to ci^o-) This equation 
implies that = a'^GMiy^jn for a' ^ Og, which is now consistent with the 
Tully-Fisher law. 

Although Milgrom's proposed modifications of gravity are consistent with 
a large amount of data, they are entirely ad hoc. Also, not only do they not 
predict the gravitational lensing by galaxies and clusters that is observed, 
it has recently been shown (Bekenstein & Sandars 1994) that in any theory 
that describes gravity by the usual tensor field of GR plus one or more scalar 
fields, the bending of light cannot exceed that predicted by GR just including 
the actual matter. Thus, if Milgrom's nonrelativistic theory, in which by 
assumption there is no dark matter, were extended to any scalar-tensor 
gravity theory, the light bending could only be due to the visible mass. 
However, the evidence is becoming increasingly convincing that the mass 
indicated by gravitational lensing in clusters of galaxies is at least as large 
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as that implied by the velocities of the galaxies and the temperature of the 
gas in the clusters (sec e.g. Wu & Fang 1996). 

Moreover, it is difficult to fit an force law into the larger framework of 
either cosmology or theoretical physics. The cosmological difficulty is that 
an force never saturates: distant masses are more important than nearby 
masses. Regarding theoretical physics, all one needs to assume in order to 
get the weak-field limit of general relativity is that gravitation is carried by a 
massless spin- two particle (the graviton): masslessness implies the standard 
force, and then spin two implies coupling to the energy-momentum 
tensor (Weinberg 1965). In the absence of an intrinsically attractive and 
plausible theory of gravity which leads to a force law at large distances, 
it seems to be preferable by far to assume GR and take dark matter seriously, 
as done below. But until the nature of the dark matter is determined — 
e.g., by discovering dark matter particles in laboratory experiments — it is 
good to remember that there may be alternative explanations for the data. 



1.3 Age, Expansion Rate, and Cosmological Constant 

1.3.1 Age of the Universe to 

The strongest lower limits for to come from studies of the stellar populations 
of globular clusters (GCs) . Standard estimates of the ages of the oldest GCs 
are tec ~ 15 — 16 Gyr (Bolte & Hogan 1995; VandenBerg, Bolte, &; Stetson 
1996; Chaboyer ct al. 1996). A frequently quoted lower limit on the age 
of GCs is 12 Gyr (Chaboyer et al. 1996), which is then an even more 
conservative lower limit on to = tac + ^^GCj where /^tcc -^0.5 Gyr is the 
time from the Big Bang until GC formation. The main uncertainty in the 
GC age estimates comes from the uncertain distance to the GCs: a 0.25 
magnitude error in the distance modulus translates to a 22% error in the 
derived cluster age (Chaboyer 1995). (We will come back to this in the next 
paragraph.) All the other obvious ways to lower the calculated tec have 
been considered and found to have limited effects, and many non-obvious 
ideas have also been explored (VandenBerg et al. 1996). For example, 
stellar mass loss is a way of lowering tec (Willson, Bowen, & Struck-Marcell 
1987), but observations constrain the reduction in to to be less than ~ 1 Gyr 
(Shi 1995, Swenson 1995). Helium sedimentation during the main sequence 
lifetime can reduce stellar ages by ~ 1 Gyr (Chaboyer & Kim 1995, D Antona 
et al. 1997). Note that the higher primordial ^He abundance implied by the 
new Tytler et al. (1996) D/H lowers the central value of the GC ages by 
perhaps 0.5 Gyr. The usual conclusion has been that to ~ 12 Gyr is probably 
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the lowest plausible value for to, obtained by pushing many but not all the 
parameters to their limits. 

However, in spring of 1997, analyses of data from the Hipparcos 
astrometric satellite have indicated that the distances to GCs assumed in 
obtaining the ages just discussed were systematically underestimated. If this 
is true, it follows that their stars at the main sequence turnoff are brighter 
and therefore younger. Indeed, there are indications that this correction 
will be largest for the lowest-metallicity clusters that had the oldest ages 
according to the standard analysis, according to Reid (1977). His analysis, 
using a sample including 15 metal-poor stars with parallaxes determined 
to better than 12% accuracy to redefine the subdwarf main sequence, gives 
distance moduli ~ 0.3 magnitudes (~ 30%) brighter than current standard 
values for his four lowest-metallicity GCs (M13, M15, M30, and M92), and 
ages {not lower limits) of ~ 12 Gyr. The shapes of the theoretical isochrones 
(Bergbusch & Vandenberg 1992) used in previous GC age estimates (e.g., 
Bolte & Hogan 1995, Sandquist et al. 1996) arc no longer acceptable fits 
to the subdwarf data with the revised distances, although the isochrones of 
D'Antona et al. (1997) give better fits to the local subdwarfs and to the 
GCs. Another analysis (Gratton et al. 1997) uses a sample including 11 
low-metallicity non-binary subdwarf stars with Hipparcos parallaxes better 
than 10% and accurate metal abundances from high-resolution spectroscopy 
to determine the absolute location of the main sequence as a function of 
metalUcity. They then derive ages for the old GCs (MIS, M68, M92, 
NGC288, NGC6752, 47 Tuc) in their GC sample of 12.11^;^ Gyr. Their 
ages are lower both because of their 0.2 mag brighter distance moduli and 
because of their better metal determinations of cluster and field stars. 

There are systematic effects that must be taken into account in the 
accurate determination of toe, including metallicity dependence and 
reddening corrections, and various physical phenomena such as stellar 
convection and helium sedimentation whose inclusion could lower ages still 
further and perhaps also bring theoretical isochrones into better agreement 
with the GC observations. Thus, there may be a period during which 
additional data is sought and theoretical models are revised before a new 
consensus emerges regarding the GC ages. But it does appear that the older 
estimates tec ~ 15 — 16 Gyr will be revised downward substantially. For 
example, in light of the new Hipparcos data, Chaboyer et al. 1997 have 
redone their Monte Carlo analysis of the effects of varying various uncertain 
parameters, and obtained tec = 11-7 i 1.4 Gyr (1<t). 

Stellar age estimates arc also relevant to another sort of argument for 
an old, low-density universe: observation of apparently old galaxies at 
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moderately high redshift (Dunlop et al. 1996). In the most extreme example 
presented so far (Spinrad et al. 1997), galaxy LDBS 53W091 at redshift 
z = 1.55 has a rest-frame spectrum very similar to that of an F6 star, and 
the claimed minimum age of 3.5 Gyr is based on standard stellar evolution 
models and assumptions about stellar populations, reddening, etc. The 
authors point out that for 3.5 Gyr to have elapsed at z = 1.55 requires 
h < 0.45 for = 1. (Note that the constraint on h is sensitive to the 
claimed age of the galaxy. From Figure LI (a), the age of an Einstein-de 
Sitter universe at z = 1.55 is 1.60/i^^ Gyr, so for 3.0 Gyr to have elapsed by 
z = 1.55 in this cosmology imposes the less restrictive requirement h < 0.53, 
for example.) Observations of old galaxies at high redshift will certainly 
constrain cosmological parameters, especially if the assumptions that go 
into the analysis can be independently verified. However, in this case an 
independent analysis (Bruzual & Magris 1997) of the same data gives a 
much younger age of 1 to 2 Gyr for LDBS 53W091 (an age of 2 Gyr at 
z = 1.55 poses no problem for an Q = 1 cosmology as long as h < 0.8), 
which these authors regard as more reliable since, unlike the earlier authors, 
they can explain all the spectral and color data. 

Stellar age estimates are of course based on standard stellar evolution 
calculations. But the solar neutrino problem reminds us that we are not 
really sure that we understand how even our nearest star operates; and the 
sun plays an important role in calibrating stellar evolution, since it is the 
only star whose age we know independently (from radioactive dating of early 
solar system material). An important check on stellar ages can come from 
observations of white dwarfs in globular and open clusters (cf. Richer et 
al. 1995). And the two detached eclipsing binaries at the main sequence 
turn-off point recently discovered in Omega Centauri can be used both to 
measure the distance to this globular cluster accurately, and to determine 
their ages using the mass- luminosity relation (Paczynski 1996). 

What if the GC age estimates are wrong for some unknown reason? The 
only other non-cosmological estimates of the age of the universe come from 
nuclear cosmochronometry — radioactive decay and chemical evolution of 
the Galaxy — and white dwarf cooling. Cosmochronometry age estimates 
are sensitive to a number of uncertain issues such as the formation history of 
the disk and its stars, and possible actinide destruction in stars (Malaney, 
Mathews, & Dearborn 1989; Mathews &: Schramm 1993). However, an 
independent cosmochronometry age estimate of 13.8 ib 3.7 Gyr has been 
obtained for a single ultra-low-metallicity star ([Fe/H]=-3.1), based on the 
measured depletion of thorium (whose half-life is 14.2 Gyr) compared to 
stable heavy r-process elements (Cowan et al. 1997; cf. Bolte 1997, Sneden 
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et al. 1996). This method wih become very important if it is possible to 
obtain accurate measurements of r-process elements for a number of very 
low metallicity stars, and the resulting age estimates are consistent. 

Independent age estimates come from the cooling of white dwarfs in the 
neighborhood of the sun. The key observation is that there is a lower limit to 
the luminosity, and therefore also the temperature, of nearby white dwarfs; 
although dimmer ones could have been seen, none have been found. The 
only plausible explanation is that the white dwarfs have not had sufficient 
time to cool to lower temperatures, which initially led to an estimate of 
9.3 ± 2 Gyr for the age of the Galactic disk (Winget et al. 1987). Since 
there was evidence (based on the pre-Hipparcos GC distances) that the 
stellar disk of our Galaxy is about 2 Gyr younger than the oldest GCs (e.g.. 
Stetson, VandenBerg, & Bolte 1996), this in turn gave an estimate of the 
age of the universe of to ~ 11 ±2 Gyr. More recent analyses (cf. Wood 1992, 
Hernanz et al. 1994) conclude that sensitivity to disk star formation history, 
and to effects on the white dwarf cooling rates due to C/0 separation at 
crystallization and possible presence of trace elements such as ^^Ne, allow 
a rather wide range of ages for the disk of about 10 it 4 Gyr. The latest 
determination of the white dwarf luminosity function, using white dwarfs in 
proper motion binaries, leads to a somewhat lower minimum luminosity and 
therefore a somewhat higher estimate of the age of the disk of ~ 10.5]ti'5 
Gyr (Oswalt et al. 1996; cf. Chabrier 1997). 

Suppose that the old GC stellar age estimates that to ~ 13 Gyr are right, 
as we will assume in much of the rest of this chapter. Figure 1.2 shows that 
to > 13 Gyr implies that h < 0.50 for 17 = 1, and that h < 0.73 even for 
r^o cis small as 0.3 in flat cosmologies (i.e., with Qq + i}\ = 1). However, 
in view of the preliminary analyses using the new Hipparcos parallaxes and 
other new data that give strikingly lower age estimates for the oldest GCs, 
we should bear in mind that to might actually be as low as ~ 11 Gyr, which 
would allow h as high as 0.6 for 0, = 1. 



1.3.2 Hubble Parameter Hq 

The Hubble parameter Hq = lOO/i km s^^ Mpc~^ remains uncertain, 
although by less than the traditional factor of two. de Vaucouleurs long 
contended that /i ~ 1. Sandage has long contended that h ~ 0.5, 
and he and Tammann still conclude that the latest data are consistent 
with h = 0.55 ± 0.05 (Sandage 1995; Sandage & Tammann 1995, 1996; 
Tammann & Federspiel 1996). A majority of observers currently favor a 
value intermediate between these two extremes, and the range of recent 
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Hubble Parameter 
(km/s/Mpc) 

Fig. 1.2. Age of the universe to as a funetion of Hubble parameter Hq in inflation 
inspired models with fio + = 1, for several values of the present-epoch 
cosmological density parameter Oq- 

determinations has been shrinking (Kennicutt, Freedman, & Mould 1995; 
Tammann 1996; Freedman 1996). 

The Hubble parameter has been measured in two basic ways: (1) 
Measuring the distance to some nearby galaxies, typically by measuring the 
periods and luminosities of Cepheid variables in them; and then using these 
"calibrator galaxies" to set the zero point in any of the several methods 
of measuring the relative distances to galaxies. (2) Using fundamental 
physics to measure the distance to some distant object(s) directly, thereby 
avoiding at least some of the uncertainties of the cosmic distance ladder 
(Rowan-Robinson 1985). The difficulty with method (1) was that there was 
only a handful of calibrator galaxies close enough for Cepheids to be resolved 
in them. However, the success of the HST Cepheid measurement of the 
distance to MlOO (Freedman et al. 1994, Ferrarese et al. 1996) shows that 
the HST Key Project on the Extragalactic Distance Scale can significantly 
increase the set of calibrator galaxies — in fact, it already has done so. 
Adaptive optics from the ground may also be able to contribute to this effort, 
although the first published result of this approach (Pierce et al. 1994) is 
not entirely convincing. The difficulty with method (2) is that in every case 
studied so far, some aspect of the observed system or the underlying physics 
remains somewhat uncertain. It is nevertheless remarkable that the results 
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of several different methods of type (2) are rather similar, and indeed not 
very far from those of method (1). This gives reason to hope for convergence. 

1.3.2.1 Relative Distance Methods 

One piece of good news is that the several methods of measuring the 
relative distances to galaxies now mostly seem to be consistent with each 
other (Jacoby et al. 1992; Fukugita, Hogan, & Peebles 1993). These 
methods use either (a) "standard candles" or (b) empirical relations between 
two measurable properties of a galaxy, one distance-independent and the 
other distance-dependent, (a) The old favorite standard candle is Type la 
supernovae; a new one is the apparent maximum luminosity of planetary 
nebulae (Jacoby ct al. 1992). Sandagc ct al. (1996) and others (van den 
Bergh 1995, Branch et al. 1996, cf. Schaefer 1996) get low values oih k, 0.55 
from HST Cepheid distances to SN la host galaxies, including the seven SNe 
la with what Sandage et al. characterize as well-observed maxima that lie 
in six galaxies for which HST Cepheid distances are now available. But 
taking account of an empirical relationship between the SN la light curve 
shape and maximum luminosity leads to higher h = 0.65 =b 0.06 (Riess, 
Press, & Kirshner 1996) or ^ = 0.63 ± 0.03 (Hamuy et al. 1996), although 
Tammann & Sandage (1995) disagree that the increase in h can be so large, 
(b) The old favorite empirical relation used as a relative distance indicator 
is the Tully-Fisher relation between the rotation velocity and luminosity 
of spiral galaxies (and the related Faber- Jackson or I?„ — a relation). A 
newer one is based on the decrease in the fluctuations in elliptical galaxy 
surface brightness on a given angular scale as comparable galaxies are seen 
at greater distances (Tonry 1991); a new SBF survey gives h = 0.81 ± 0.06 
(Tonry et al. 1997). 

The "mid-term" value of the Hubble constant from the HST key project 
is /i = 0.73 ± 0.10 (Frcedman et al. 1997). This is based on the standard 
distance to the LMC of 50 kpc (corresponding to a distance modulus of 
18.50). But the preliminary results from the Hipparcos astrometric satellite 
suggest that the Cepheid distance scale must be recalibrated, and that the 
quoted distance to the LMC is too low by about 10% (Feast k. Catchpole 
1997, Feast &: Whitelock 1997). An increase in the LMC distance of about 
7% is also obtained using the preliminary Hipparcos recalibration of the zero 
point and metallicity dependence of the RR Lyrae distance scale (Gratton et 
al. 1997, Ried 1997; cf. Alcock et al. 1996b), thus removing a long-standing 
discrepancy. The implication is that the Hubble parameter determined by 
Cepheid calibrators must be decreased, by perhaps 10%. This applies to 
the HST key project, and it also applies to the SN la results for /i, which 
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are based on Cepheid distances; thus, for example, the Hamuy et al. (1996) 
value would decrease to about h = 0.57, with a corresponding to = 11.4 Gyr 

for n = i. 

1.3.2.2 Fundamental Physics Approaches 

The fundamental physics approaches involve either Type la or Type II 

supernovac, the Sunyaev-Zel'dovich (S-Z) effect, or gravitational lensing. 
All are promising, but in each case the relevant physics remains somewhat 
uncertain. 

The ^^Ni radioactivity method for determining Hq using Type la SN 
avoids the uncertainties of the distance ladder by calculating the absolute 
luminosity of Type la supernovae from first principles using plausible but as 
yet unproved physical models. The first result obtained was that h = 0.61 ± 
0.10 (Arnet, Branch, & Wheeler 1985; Branch 1992); however, another study 
(Leibundgut &; Pinto 1992; cf. Vaughn et al. 1995) found that uncertainties 
in extinction (i.e., light absorption) toward each supernova increases the 
range of allowed h. Demanding that the ^^Ni radioactivity method agree 
with an expanding photosphere approach leads to h = 0.60to (Nugent et 
al. 1995). The expanding photosphere method compares the expansion rate 
of the SN envelope measured by redshift with its size increase inferred from 
its temperature and magnitude. This approach was first applied to Type II 
SN; the 1992 result h = 0.6±0.1 (Schmidt, Kirschner, & Eastman 1992) was 
subsequently revised upward by the same authors to h = 0.73 it 0.06 ± 0.07 
(1994). However, there are various complications with the physics of the 
expanding envelope (Ruiz-Lapuente et al. 1995; Eastman, Schmidt, & 
Kirshner 1996). 

The S-Z effect is the Compton scattering of microwave background 
photons from the hot electrons in a foreground galaxy cluster. This can 

be used to measure Hq since properties of the cluster gas measured via the 
S-Z effect and from X-ray observations have different dependences on Hq. 
The result from the first cluster for which sufficiently detailed data was 
available, A665 (at z = 0.182), was h = (0.4 - 0.5) ± 0.12 (Birkinshaw, 
Hughes, Sz Arnoud 1991); combining this with data on A2218 {z = 0.171) 
raised this somewhat to h = 0.55 it 0.17 (Birkinshaw & Hughes 1994). Early 
results from the ASCA X-ray satellite gave h = 0.47 ± 0.17 for A665 and 
h = 0.4llgj| for CL0016+16 {z = 0.545) (Yamashita 1994). A few S-Z 
results have been obtained using millimeter-wave observations (Wilbanks 
1994), and this method may allow more such measurements soon. New 
results for A2218 and A1413 (z = 0.14) using the Ryle radio telescope 
and ROSAT X-ray data gave h = O.SStoH and h = O.ArtoH, respectively 
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(Lasenby 1996). New results from the OVRO 5.5m telescope for the four 
X-ray brightest clusters give h = 0.54ib0.14 (Myers ct al. 1997). Corrections 
for the near-relativistic electron motions (Rcphacli 1995) and for Icnsing by 
the cluster (Loeb &: Refregier 1996) may raise these estimates for Hq a 
little, but it seems clear that the S-Z results favor a smaller value than 
many optical astronomers obtain. However, since the S-Z measurement of 
Hq is affected by the isothcrmality of the clusters (Roettiger et al. 1997) 
and the unknown orientation of the chistcr cllipticity with respect to the 
line of sight, and the errors in the derived values remain rather large, this 
lower S-Z Hq can only become convincing with more detailed observations 
and analyses of a significant number of additional clusters. Perhaps this will 
be possible within the next several years. 

Several quasars have been observed to have multiple images separated by a 
few arc seconds; this phenomenon is interpreted as arising from gravitational 
lensing of the source quasar by a galaxy along the line of sight. In the 
first such system discovered, QSO 0957+561 (z = 1.41), the time delay 
At between arrival at the earth of variations in the quasar's luminosity 
in the two images has been measured to be, e.g., 409 =b 23 days (Pelt et 
al. 1994), although other authors found a value of 540 it 12 days (Press, 
Rybicki, & Hewitt 1992). The shorter At has now been confirmed by 
the observation of a sharp drop in Image A of about 0.1 mag in late 
December 1994 (Kundic et al. 1995) followed by a similar drop in Image 
B about 405-420 days later (Kundic et al. 1997a). Since At ^ O'^H^^, 
this observation allows an estimate of the Hubble parameter, with the early 
results h = 0.50 ± 0.17 (Rhee 1991), ot h = 0.63 ± 0.21 {h = 0.42 ± 0.14) 
including (neglecting) dark matter in the lensing galaxy (Roberts et al. 
1991), with additional uncertainties associated with possible microlensing 
and unknown matter distribution in the lensing galaxy and the cluster 
in which this is the first-ranked galaxy. Deep images allowed mapping 
of the gravitational potential of the cluster (at z = 0.36) using weak 
gravitational lensing, which led to the conclusion that h < 0.70(1. lyr/At) 
(Dahle, Maddox, & Lilje 1994; Rhee et al. 1996, Fischer et al. 1997). 
Detailed study of the lenscd QSO images (which include a jet) constrains the 
lensing and implies h = 0.85(1 — K)(l.lyr/At) < 0.85, where the upper limit 
follows because the convergence due to the cluster At > 0, or alternatively 
h = 0.85(c7/322kms~^)2(l.lyr/At) without uncertainty concerning the 
cluster if the one-dimensional velocity dispersion a in the core of the giant 
elliptical galaxy responsible for the lensing can be measured (Grogin & 
Narayan 1996). The latest results for h from 0957+561, using all available 
data, are h = 0.64 + 0.13 (95% C.L.) (Kundic et al. 1997a), h = 0.62 + 0.07 
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(Falco et al. 1997, where the error does not include systematic errors in the 
assumed form of the mass distribution in the lens; uncertainties can also be 
reduced with new HST images of the system, allowing improved accuracy 
in the lens galaxy position). 

The first quadruple- image quasar system discovered was PGl 115+080. 
Using a recent series of observations (Schechter et al. 1997), the time delay 
between images B and C has been determined to be about 24 it 3 days, 
or 25l3'g days by an alternative analysis (BarKana 1997). A simple model 
for the Icnsing galaxy and the nearby galaxies then leads to h = 0.42 ± 
0.06 (Schechter et al. 1997) oi h = 0.41 ± 0.12 (95% C.L.) (BarKana, 
private communication), although higher values for h are obtained by a 
more sophisticated analysis: h = 0.60 it 0.17 (Keeton & Kochanek 1996), 
h = 0.52±0.14 (Kundic et al. 1997b). The results depend on how the lensing 
galaxy and those in the compact group of which it is a part are modelled. 
Such models need to be constrained by new HST observations, especially of 
the light profile in the lensing galaxy, and spectroscopy to better determine 
the velocity dispersion of the lensing galaxy and of the group. 

Although the most recent time-delay results for h from both lensed 
quasar systems are remarkably close, the uncertainty in the h determination 
by this method remains rather large. But it is reassuring that this 
completely independent method gives results consistent with the other 
determinations. The time-delay method is promising (Blandford &: Kundic 
1996), and when these systems are better understood and/or delays are 
reliably measured in several other multiple-image quasar systems, such as 
B1422+231 (Hammer, Rigaut, & Angonin-Willaimc 1995, Hjorth et al. 
1996), or radio Einstein-ring systems, such as PKS 1830-211 (van Ommen 
et al. 1995) or B0218-I-357 (Corbett et al. 1996), that should lead to a more 
precise and reliable value for Hq. 

1.3.2.3 Correcting for Virgocentric Infall 

What about the HST Cepheid measurement of Hq, giving h = 0.80 ± 0.17 
(Preedman et al. 1994), which received so much attention in the press? This 
calculated value is based on neither of the two methods (A) or (B) above, 
and it should not be regarded as being very reliable. Instead this result is 
obtained by assuming that MlOO is at the core of the Virgo cluster, and 
dividing the sum of the recession velocity of Virgo, about 1100 km s~^, plus 
the calculated "infall velocity" of the local group toward Virgo, about 300 
km s~^, by the measured distance to MlOO of 17.1 Mpc. (These recession and 
infall velocities are both a little on the high side, compared to other values 
one finds in the literature.) Adding the "infall velocity" is necessary in this 
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method in order to correct the Virgo recession velocity to what it would 
be were it not for the gravitational attraction of Virgo for the Local Group 
of galaxies, but the problem with this is that the net motion of the Local 
Group with respect to Virgo is undoubtedly affected by much besides the 
Virgo cluster — e.g., the "Great Attractor." For example, in our CHDM 
supercomputer simulations (which appear to be a rather realistic match 
to observations), galaxies and groups at about 20 Mpc from a Virgo-sized 
cluster often have net outflowing rather than infalling velocities. Note that 
if the net "infall" of MlOO were smaller, or if MlOO were in the foreground 
of the Virgo cluster (in which case the actual distance to Virgo would be 
larger than 17.1 Mpc), then the indicated Hq would be smaller. 

Freedman et al. (1994) gave an alternative argument that avoids the 
"infall velocity" uncertainty: the relative galaxy luminosities indicate that 
the Coma cluster is about six times farther away than the Virgo cluster, and 
peculiar motions of the Local Group and the Coma cluster are relatively 
small corrections to the much larger recession velocity of Coma; dividing 
the recession velocity of the Coma cluster by six times the distance to MlOO 
again gives Hq « 80. However, this approach still assumes that MlOO is in 
the core rather than the foreground of the Virgo cluster; and in deducing the 
relative distance of the Coma and Virgo clusters it assumes that the galaxy 
luminosity functions in each are comparable, which is uncertain in view 
of the very different environments. More general arguments by the same 
authors (Mould et al. 1995) lead them to conclude that h = 0.73 ± 0.11 
regardless of where MlOO lies in the Virgo cluster. But Tammann et al. 
(1996), using all the available HST Cepheid distances and their own complete 
sample of Virgo spirals, conclude that h « 0.54. 

1.3.2.4 Conclusions on Hq 

To summarize, many observers, using mainly relative distance methods, 
favor a value h ~ 0.6 — 0.8 although Sandage's group and some others 
continue to get h « 0.5 — 0.6 and all of these values may need to be 
reduced by something like 10% if the full Hipparcos data set bears out the 
preliminary reports discussed above. Meanwhile the fundamental physics 
methods typically lead to /i ~ 0.4 — 0.7. Among fundamental physics 
approaches, there has been important recent progress in measuring h via 
time delays between different images of gravitationally lensed quasars, with 
the latest analyses of both of the systems with measured time delays giving 
0.6 ±0.1. 

The fact that the fundamental physics measurements giving lower 
values for h (via time delays in gravitationally lensed quasars and the 
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Sunyaev-Zel'dovich effect) are mostly of more distant objects lias suggested 
to some authors (Turner, Cen, Sz Ostriker 1992; Wu et al. 1996) that 
the local universe may actually be underdense and therefore be expanding 
faster than is typical. But in reasonable models where structure forms from 
Gaussian fluctuations via gravitational instability, it is extremely unlikely 
that a sufficiently large region has a density sufficiently smaller than average 
to make more than a rather small difference in the value of h measured 
locally (Suto, Suginohara, &; Inagaki 1995; Shi &: Turner 1997). Moreover, 
the small dispersion in the corrected maximum luminosity of distant Type 
la supernovae found by the LBL Supernova Cosmology Project (Kim et 
al. 1997) compared to nearby SNe la shows directly that the local and 
cosmological values of Hq are approximately equal. The maximum deviation 
permitted is about 10%. Interestingly, preliminary results using 44 nearby 
Type la supernovae as yardsticks suggest that the actual deviation is about 
5-7%, in the sense that in our local region of the universe, out to a radius 
of about 70/i~^Mpc (the distance of the Northern Great Wall), Hq is this 
much larger than average (A. Dekel, private communication). The combined 
effect of this and the Hipparcos correction would, for example, reduce the 
"mid-term" value h ~ 0.73 from the EST Key Project on the Extragalactic 
Distance Scale, to /i ~ 0.63. 

There has been recent observational progress in both relative distance and 
fundamental physics methods, and it is likely that the Hubble parameter will 
be known reliably to 10% within a few years. Most recent measurements are 
consistent with /i = 0.6 it 0.1, corresponding to a range to = 6.52/i~^Gyr = 
9.3—13.0 Gyr for 0, = 1 — in good agreement with the preliminary estimates 
of the ages of the oldest globular clusters based on the new data from the 
Hipparcos astrometric satellite. 

1.3.3 Cosmological Constant A 

Inflation is the only known solution to the horizon and flatness problems and 
the avoidance of too many GUT monopoles. And inflation has the added 
bonus that at no extra charge (except the perhaps implausibly fine-tuned 
adjustment of the self-coupling of the inflaton field to be adequately 
small), simple inflationary models predict a near-Zel'dovich primordial 
spectrum (i.e., Pp{k) oc k'^" with rip « 1) of adiabatic Gaussian primordial 
fluctuations — which seems to be consistent with observations. All simple 
inflationary models predict that the curvature is vanishingly small, although 
inflationary models that are extremely contrived (at least, to my mind) 
can be constructed with negative curvature and therefore ~ 1 without a 
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cosmological constant (see § 1.6. 6| below). Thus most authors who consider 



inflationary models impose the condition A; = 0, or JIq + = 1 where 
= A/{3Hq). This is what is assumed in ACDM models, and it is 
what was assumed in Fig. |1.2| . (Note that is used to refer only to the 
density of matter and energy, not including the cosmological constant, whose 
contribution in 0, units is 

The idea of a nonvanishing A is commonly considered unattractive. There 
is no known physical reason why A should be so small {Q\ = 1 corresponds 
to PA ~ 10~^^ eV*, which is small from the viewpoint of particle physics), 
though there is also no known reason why it should vanish (cf. Weinberg 
1989, 1996). A very unattractive feature of A 7^ cosmologies is the fact 
that A must become important only at relatively low redshift — why not 
much earlier or much later? Also Oa ~ ^^o implies that the universe has 
recently entered an inflationary epoch (with a de Sitter horizon comparable 
to the present horizon). The main motivations for A > cosmologies are (1) 
reconciling inflation with observations that seem to imply < li and (2) 
avoiding a contradiction between the lower limit to ~ 13 Gyr from globular 
clusters and to = {2/3)Hq^ = 6.52/i~^ Gyr for the standard 17 = 1, A = 
Einstein-de Sitter cosmology, if it is really true that h > 0.5. 

The cosmological effects of a cosmological constant are not difficult to 
understand (Lahav et al. 1991; Carroll, Press, &: Turner 1992). In the 
early universe, the density of energy and matter is far more important than 
the A term on the r.h.s. of the Friedmann equation. But the average 
matter density decreases as the universe expands, and at a rather low 
redshift {z ~ 0.2 for f^o = 0.3) the A term finally becomes dominant. If 
it has been adjusted just right, A can almost balance the attraction of 
the matter, and the expansion nearly stops: for a long time, the scale 
factor a = {1 + z)~^ increases very slowly, although it ultimately starts 
increasing exponentially as the universe starts inflating under the influence 
of the increasingly dominant A term (see Fig. LI). The existence of a 
period during which expansion slows while the clock runs explains why to 
can be greater than for A = 0, but this also shows that there is an increased 
likelihood of finding galaxies at the redshift interval when the expansion 
slowed, and a correspondingly increased opportunity for lensing of quasars 
(which mostly lie at higher redshift z ^ 2) by these galaxies. 

The frequency of such lensed quasars is about what would be expected 
in a standard = 1, A = cosmology, so this data sets fairly stringent 
upper limits: < 0.70 at 90% C.L. (Maoz & Rix 1993, Kochanek 1993), 
with more recent data giving even tighter constraints: Q\ < 0.66 at 95% 
confidence if f^o + = 1 (Kochanek 1996b). This limit could perhaps be 
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weakened if there were (a) significant extinction by dust in the E/SO galaxies 
responsible for the lensing or (b) rapid evolution of these galaxies, but there 
is much evidence that these galaxies have little dust and have evolved only 
passively for z ^ 1 (Steidel, Dickinson, &: Persson 1994; Lilly et al. 1995; 
Schade et al. 1996). (An alternative analysis by Im, Griffiths, & Ratnatunga 
1997 of some of the same optical lensing data considered by Kochanek 1996b 
leads them to deduce a value fl^ = O.GA'^q'Iq, which is barely consistent with 
Kochanek's upper limit. A recent paper — Malhotra, Rhodes, &: Turner 
1997 — presents evidence for extinction of quasars by foreground galaxies 
and claims that this weakens the lensing bound to 0,a < 0.9, but there is 
no quantitative discussion in the paper to justify this claim. Mailer, Flores, 
& Primack 1997 shows that edge-on disk galaxies can lens quasars very 
effectively, and discusses a case in which optical extinction is significant. 
But the radio observations discussed by Falco, Kochanek, & Munoz 1997, 
which give a 2a limit Q\ < 0.73, will not be affected by extinction.) 

Yet another constraint comes from number counts of bright E/SO galaxies 
in HST images (Driver et al. 1996), since as was just mentioned these 
galaxies appear to have evolved rather little since z ~ 1. The number counts 
are just as expected in the J7 = 1, A = Einstein-de Sitter cosmology. Even 
allowing for uncertainties due to evolution and merging of these galaxies, this 
data would allow Q\ as large as 0.8 in flat cosmologies only in the unlikely 
event that half the Sa galaxies in the deep HST images were misclassified as 
E/SO. This number-count approach may be very promising for the future, 
as the available deep HST image data and our understanding of galaxy 
evolution both increase. 

A model-dependent constraint comes from a detailed simulation of 
ACDM (Klypin, Primack, & Holtzman 1996, hereafter KPH96): a 
COBE-normalized model with 0,q = 0.3, JIa = 0.7, and h = 0.7 has far too 
much power on small scales to be consistent with observations, unless there 
is unexpectedly strong scale-dependent antibiasing of galaxies with respect 
to dark matter. (This is discussed in more detail in § 1.7.4 below.) For 
ACDM models, the simplest solution appears to be raising Qq, lowering Hq, 
and tilting the spectrum (up < 1), though of course one could alternatively 
modify the primordial power spectrum in other ways. 

Finally, from a study of their first seven high-redshift Type la supernovae, 
Perlmutter et al. (1996) deduced that Q\ < 0.51 at the 95% confidence level. 



(This is discussed in more detail in § 1.4.1, just below. 



Figure 1.2 shows that with ^ 0.7, the cosmological constant does not 
lead to a very large increase in Iq compared to the Einstein-de Sitter case, 
although it may still be enough to be significant. For example, the constraint 
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that Iq > 13 Gyr requires h < 0.5 for = 1 and A = 0, but this becomes 
h < 0.70 for flat cosmologies with < 0.66. 



1.4 Measuring $7o 

The present author, hke many theorists, regards the Einstein-de Sitter (fi = 
1, A = 0) cosmology as the most attractive one. For one thing, there are 
only three possible constant values for O — 0, 1, and oo — of which the 
only one that can describe our universe is = 1. Also, as will be discussed 



in more detail in § 1.6.2, cosmic inflation is the only known solution for 



several otherwise intractable problems, and all simple inflationary models 
predict that the universe is flat, i.e. that + ^^a = 1- Since there is 
no known physical reason for a non-zero cosmological constant, and as just 
discussed in § 1.3.3| there are strong observational upper limits on it (e.g.. 



from gravitational lensing), it is often said that inflation favors 17 = 1. 

Of course, theoretical prejudice is not necessarily a reliable guide. In 
recent years, many cosmologists have favored rig ~ 0.3, both because of 
the Hq — to constraints and because cluster and other relatively small-scale 
measurements have given low values for Qq. (For a recent summary of 
arguments favoring low ilo ^ 0.2 and A = 0, see Coles & Ellis 1997; see also 
the chapters by Bahcall and Peebles in this book. A recent review which 
notes that larger scale measurements favor higher 0,q is Dekel, Burstein, &; 
White 1997.) However, in light of the new Hipparcos data, the Hq — to 
data no longer so strongly disfavor $7 = 1. Moreover, as is discussed in 
more detail below, the small-scale measurements are best regarded as lower 
limits on Qq- At present, the data does not permit a clear decision whether 
fio ~ 0.3 or 1, but there are promising techniques that may give definitive 
measurements soon. 



1.4 -1 Very Large Scale Measurements 

Although it would be desirable to measure find A through their effects on 
the large-scale geometry of space-time, this has proved difficult in practice 
since it requires comparing objects at higher and lower redshift, and it is 
hard to separate selection effects or the effects of the evolution of the objects 
from those of the evolution of the universe. For example, Kellermann (1993), 
using the angular-size vs. redshift relation for compact radio galaxies, 
obtained evidence favoring O ~ 1; however, selection effects may invalidate 
this approach (Dabrowski, Lasenby, & Saunders 1995). To cite another 
example, in "redshift-volume" tests (e.g. Loh & Spillar 1986) involving 
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number counts of galaxies per redshift interval, how can we tell whether the 
galaxies at redshift z ~ 1 correspond to those at z ~ 0? Several galaxies 
at higher redshift might have merged, and galaxies might have formed or 
changed luminosity at lower redshift. Eventually, with extensive surveys of 
galaxy properties as a function of redshift using the largest telescopes such 
as Keck, it should be possible to perform classical cosmological tests at least 
on particular classes of galaxies — that is one of the goals of the Keck DEEP 
project. Geometric effects are also on the verge of detection in small-angle 



cosmic microwave background (CMB) anisotropies (see § 1.4.8 and § 1.7.4| ). 

At present, perhaps the most promising technique involves searching 
for Type la supernovae (SNe la) at high-redshift, since these are the 
brightest supernovae and the spread in their intrinsic brightness appears to 
be relatively small. Perlmutter et al. (1996) have recently demonstrated 
the feasibility of finding significant numbers of such supernovae, but a 
dedicated campaign of follow-up observations of each one is required in 
order to measure Q,q by determining how the apparent brightness of the 
supernovae depends on their redshift. This is therefore a demanding 
project. It initially appeared that ~ 100 high redshift SNe la would be 
required to achieve a 10% measurement of = — However, 

using the correlation mentioned earlier between the absolute luminosity of 
a SN la and the shape of its light curve (slower decline correlates with 
higher peak luminosity), it now appears possible to reduce the number 
of SN la required. The Perlmutter group has now analyzed seven high 
redshift SN la by this method, with the result for a flat universe that 
Q.0 = I - 0.^ = 0.94+[j:^|, or equivalently 0.^ = 0.06l:[j:^^ (< 0.51 at the 
95% confidence level) (Perlmutter et al. 1996). For a A = cosmology, 
they find = O.SSlg gg. In November 1995 they discovered an additional 
11 high-redshift SN la, and they have subsequently discovered many more. 
Other groups, collaborations from ESO and MSSSO/CfA/CTIO, are also 
searching successfully for high-redshift supernovae to measure VLq (Garnavich 
et al. 1996). There has also been recent progress understanding the physical 
origin of the SN la luminosity-light curve correlation, and in discovering 
other such correlations. At the present rate of progress, a reliable answer 
may be available within perhaps a year or two if a consensus emerges from 
these efforts. 



1.4-2 Large-scale Measurements 

f^o has been measured with some precision on a scale of about ~ 50 /i~^Mpc, 
using the data on peculiar velocities of galaxies, and on a somewhat larger 
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scale using redshift surveys based on the IRAS galaxy catalog. Since the 
results of all such measurements to date have been reviewed in detail (see 
Dekel 1994, Strauss & Willick 1995, and Dekel's chapter in this volume), only 
brief comments arc provided here. The "POTENT" analysis tries to recover 
the scalar velocity potential from the galaxy peculiar velocities. It looks 
reliable, since it reproduces the observed large scale distribution of galaxies 

— that is, many galaxies are found where the converging velocities indicate 
that there is a lot of matter, and there are voids in the galaxy distribution 
where the diverging velocities indicate that the density is lower than average. 
The comparison of the IRAS redshift surveys with POTENT and related 
analyses typically give fairly large values for the parameter /?/ = D,Q'^/bi 
(where bj is the biasing parameter for IRAS galaxies), corresponding to 
0.3 ^ ^ 3 (for an assumed bj = 1.15). It is not clear whether it will be 
possible to reduce the spread in these values significantly in the near future 

— probably both additional data and a better understanding of systematic 
and statistical effects will be required. 

A particularly simple way to deduce a lower limit on Qq from the 
POTENT peculiar velocity data was proposed by Dekel & Rees (1994), 
based on the fact that high-velocity outflows from voids are not expected in 
low-O models. Data on just one void indicates that 0,q > 0.3 at the 97% C.L. 
This argument is independent of assumptions about A or galaxy formation, 
but of course it does depend on the success of POTENT in recovering the 
peculiar velocities of galaxies. 

However, for the particular cosmological models that are at the focus 
of this review — CHDM and ACDM — stronger constraints are available. 
This is because these models, in common with almost all CDM variants, 
assume that the probability distribution function (PDF) of the primordial 
fluctuations was Gaussian. Evolution from a Gaussian initial PDF to 
the non-Gaussian mass distribution observed today requires considerable 
gravitational nonlinearity, i.e. large 0,. The PDF deduced by POTENT from 
observed velocities (i.e., the PDF of the mass, if the POTENT reconstruction 
is reliable) is far from Gaussian today, with a long positive-fluctuation tail. 
It agrees with a Gaussian initial PDF if and only if is about unity or 
larger: f^o < 1 is rejected at the 2a level, and 0,q < 0.3 is ruled out at > 4o" 
(Nusser & Dekel 1993; cf. Bernardeau et al. 1995). 



1.4-3 Measurements on Scales of a Few Mpc 

On smaller length scales, there are many measurements that are 
consistent with a smaller value of (e-g- Peebles 1993, esp. §20). 
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For example, the cosmic virial theorem gives lh~^ Mpc) Ri 

0.15[(t(1/i~^ Mpc)/(300kms~"^)]^, where a{lh~^Mpc) here represents the 
relative velocity dispersion of galaxy pairs at a separation of Ih^^ Mpc. 
Although the classic paper (Davis &: Peebles 1983) which first measured 
a{lh~^ Mpc) using a large redshift survey (CfAl) got a value of 340 km s~^, 
this result is now known to be in error since the entire core of the Virgo 
cluster was inadvertently omitted (Somerville, Davis, & Primack 1996); if 
Virgo is included, the result is ~ 500 - 600 km s"^ (cf. Mo et al. 1993, 
Zurek et al. 1994), corresponding to 0(~ lh~^ Mpc) « 0.4 — 0.6. Various 
redshift surveys give a wide range of values for a{lh-^ Mpc) ~ 300 - 750 
km s~^, with the most salient feature being the presence or absence of rich 
clusters of galaxies; for example, the IRAS galaxies, which are not found 
in clusters, have a{lh~^Mpc) « 320 km s~^ (Fisher et al. 1994), while 
the northern CfA2 sample, with several rich clusters, has much larger a 
than the SSRS2 sample, with only a few relatively poor clusters (Marzke et 
al. 1995; Somerville, Primack, & Nolthenius 1996). It is evident that the 
a{lh~^ Mpc) statistic is not a very robust one. Moreover, the finite sizes of 
the dark matter halos of galaxies and groups complicates the measurement 
of n using the CVT, generally resulting in a significant underestimate of the 
actual value (Bartlctt &; Blanchard 1996, Suto & Jing 1996). 

A standard method for estimating 0, on scales of a few Mpc is based 
on applying virial estimates to groups and clusters of galaxies to try to 
deduce the total mass of the galaxies including their dark matter halos from 
the velocities and radii of the groups; roughly, GM ~ rv^. (What one 
actually does is to pretend that all galaxies have the same mass-to-light 
ratio M/L, given by the median M/L of the groups, and integrate over the 
luminosity function to get the mass density (Kirschner, Oemler, &; Schechter 
1979; Huchra & Geller 1982; Ramella, Geller, & Huchra 1989). The typical 
result is that ri(~ lh~^Mpc) ~ 0.1 — 0.2. However, such estimates are at 
best lower limits, since they can only include the mass within the region 
where the galaxies in each group can act as test particles. It has been 
found in CHDM simulations (Nolthenius, Klypin, & Primack 1997) that the 
effective radius of the dark matter distribution associated with galaxy groups 
is typically 2-3 times larger than that of the galaxy distribution. Moreover, 
we find a velocity biasing (Carlberg & Couchman 1989) factor in CHDM 
groups b':!,''^ = Wgai,rms/^^DM,rms ~ 0.75, whose inverse squared enters in the Q 
estimate. Finally, we find that groups and clusters are typically elongated, 
so only part of the mass is included in spherical estimators. These factors 
explain how it can be that our = 1 CHDM simulations produce group 
velocity dispersions that are fully consistent with those of observed groups. 



32 



Dark Matter and Structure Formation 



even with statistical tests such as the median rms internal group velocity 
vs. the fraction of galaxies grouped (Nolthenius, Klypin, & Primack 1994, 
1997). This emphasizes the point that local estimates of Q, are at best lower 
limits on its true value. 

However, a new study by the Canadian Network for Observational 
Cosmology (CNOC) of 16 clusters at z ~ 0.3 mostly chosen from the 
Einstein Medium Sensitivity Survey (Henry et al. 1992) was designed to 
allow a self-contained measurement of Qq from a field AI/ L which in turn was 
deduced from their measured cluster M/L. The result was fio = 0.19 it 0.06 
(Carlberg et al. 1997a,c). These data were mainly compared to standard 
CDM models, and they probably exclude Q = 1 in such models. But it 
remains to be seen whether alternatives such as a mixture of cold and hot 
dark matter could fit the data. 

Another approach to estimating Q from information on relatively small 
scales has been pioneered by Peebles (1989, 1990, 1994). It is based on 
using the least action principle (LAP) to reconstruct the trajectories of the 
Local Group galaxies, and the assumption that the mass is concentrated 
around the galaxies. This is perhaps a reasonable assumption in a low-fi 
universe, but it is not at all what must occur in an J7 = 1 universe where 
most of the mass must lie between the galaxies. Although comparison 
with Q = I N-body simulations showed that the LAP often succeeds 
in qualitatively reconstructing the trajectories, the mass is systematically 
underestimated by a large factor by the LAP method (Branchini &; Carlberg 
1994). Surprisingly, a different study (Dunn & Laflamme 1995) found 
that the LAP method underestimates by a factor of 4-5 even in an 
ilo = 0.2 simulation; the authors say that this discrepancy is due to the 
LAP neglecting the effect of "orphans" — dark matter particles that are 
not members of any halo. Shaya, Peebles, and Tully (1995) have recently 
attempted to apply the LAP to galaxies in the local supercluster, again 
getting low f^o- The LAP approach should be more reliable on this larger 
scale, but the method still must be calibrated on N-body simulations of both 
high- and Iow-JIq models before its biases can be quantified. 



1.4-4 Estimates on Galaxy Halo Scales 

A classic paper by Little &; Tremaine (1987) had argued that the available 
data on the Milky Way satellite galaxies required that the Galaxy's halo 
terminate at about 50 kpc, with a total mass of only about 2.5 x lO^^M©. 
But by 1991, new data on local satellite galaxies, especially Leo I, became 
available, and the Little- Tremaine estimator increased to 1.25 x lO^^M©. A 
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recent, detailed study finds a mass inside 50 kpc of (5.4 it 1.3) x 10^^ Mq 
(Kochanck 1996a). 

Work by Zaritsky et al. (1993) has shown that other spiral galaxies 
also have massive halos. They collected data on satellites of isolated spiral 
galaxies, and concluded that the fact that the relative velocities do not fall off 
out to a separation of at least 200 kpc shows that massive halos are the norm. 
The typical rotation velocity of 200 — 250 km s~^ implies a mass within 
200 kpc of ~ 2 X IO^^Mq. A careful analysis taking into account selection 
effects and satellite orbit uncertainties concluded that the indicated value of 
fio exceeds 0.13 at 90% confidence (Zaritsky & White 1994), with preferred 
values exceeding 0.3. Newer data suggesting that relative velocities do not 
fall off out to a separation of ~ 400 kpc (Zaritsky et al. 1997) presumably 
would raise these estimates. 

However, if galaxy dark matter halos are really so extended and massive, 
that would imply that when such galaxies collide, the resulting tidal tails of 
debris cannot be flung very far. Therefore, the observed merging galaxies 
with extended tidal tails such as NGC 4038/39 (the Antennae) and NGC 
7252 probably have halo:(disk+bulge) mass ratios less than 10:1 (Dubinski, 
Mihos, Sz Hernquist 1996), unless the stellar tails are perhaps made during 
the collision process from gas that was initially far from the central galaxies 
(J. Ostriker, private communication, 1996); the latter possibility can be 
checked by determining the ages of the stars in these tails. 

A direct way of measuring the mass and spatial extent of many galaxy 
dark matter halos is to look for the small distortions of distant galaxy images 
due to gravitational lensing by foreground galaxies. This technique was 
pioneered by Tyson et al. (1984). Though the results were inconclusive 
(Kovner &; Milgrom 1987), powerful constraints could perhaps be obtained 
from deep HST images or ground-based images with excellent seeing. Such 
fields would also be useful for measuring the correlated distortions of galaxy 
images from large-scale structure by weak gravitational lensing; although a 
pilot project (Mould et al. 1994) detected only a marginal signal, a reanalysis 
detected a significant signal suggesting that JloCs ~ 1 (Villumsen 1995). 
Several groups are planning major projects of this sort. The first results 
from an analysis of the Hubble Deep Field gave an average galaxy mass 
interior to 20/1"^ kpc of 5.9^2,^ x 10^^ /j-^Mq (Dell' Antonio & Tyson 1996). 

1.4-5 Cluster Baryons vs. Big Bang Nucleosynthesis 

A review (Copi, Schramm, &: Turner 1995) of Big Bang Nucleosynthesis 
(BBN) and observations indicating primordial abundances of the light 
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isotopes concludes that 0.009^"^ < < 0.02h~^ for concordance with 
all the abundances, and 0.006/i~^ < < 0.03h~^ if only deuterium is 
used. For h = 0.5, the corresponding upper limits on ^Ijy arc 0.08 and 
0.12, respectively. The observations (Songaila et al. 1994a, Carswell et al. 
1994) of a possible deuterium line in a hydrogen cloud at redshift z = 3.32 
in the spectrum of quasar 0014+813, indicating a deuterium abundance 
D/H~ 2 X 10~^ (and therefore fi^ < 0.006/i~^), are inconsistent with D/H 
observations by Tytlcr and collaborators (Tytlcr ct al. 1996, Buries &; Tytler 
1996) in systems at z = 3.57 (toward Q1937-1009) and at z = 2.504, but with 
a deuterium abundance about ten times lower. These lower D/H values are 
consistent with solar system measurements of D and ^Hc, and they imply 
Qf^h'^ = 0.024 ± 0.05, or in the range 0.08-0.11 for h = 0.5. If these 
represent the true D/H, then if the earlier observations were correct they 
were most probably of a Lya forest line. Rugers k, Hogan (1996) argue that 
the width of the z = 3.32 absorption features is better fit by deuterium, 
although they admit that only a statistical sample of absorbers will settle 
the issue. There is a new possible detection of D at z = 4.672 in the 
absorption spectrum of QSO BR1202-0725 (Wampler et al. 1996) and at 
z = 3.086 toward Q0420-388 (Carswell 1996), but they can only give upper 
limits on D/H. Wampler (1996) and Songaila ct al. (1997) claim that Tytler 
et al. (1996) have overestimated the HI column density in their system, 
and therefore underestimated D/H. But Buries &: Tytler (1996) argue that 
the two systems that they have analyzed are much more convincing as 
real detections of deuterium, that their HI column density measurement 
is reliable, and that the fact that they measure the same D/H~ 2.4 x 10~^ 
in both systems makes it likely that this is the primordial value. Moreover, 
Tytler, Buries, &; Kirkman (1996) have recently presented a higher resolution 
spectrum of Q0014+813 in which "deuterium absorption is neither required 
nor suggested," which would of course completely undercut the argument 
of Hogan and collaborators for high D/H. Finally, the Tytler group has 
analyzed their new Keck LRIS spectra of the absorption system toward 
Q1937-1009, and they say that the lower HI column density advocated by 
Songaila et al. (1997) is ruled out (Buries and Tytler 1997). Of course, one 
or two additional high quality D/H measurements would be very helpful to 
really settle the issue. 

There is an entirely different line of argument that also favors the higher 
implied by the lower D/H of Tytler et al. This is the requirement that 
the high-redshift intergalactic medium contain enough neutral hydrogen to 
produce the observed Lymanct forest clouds given standard estimates of the 
ultraviolet ionizing flux from quasars. The minimum required O;, ^ 0.05/i^q^ 
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(Gnedin &: Hui 1996, Weinberg et al. 1997) is considerably higher than that 
advocated by higher D/H values, but consistent with that implied by the 
lower D/H measurements. 

Yet another argument favoring the D/H of Tytler et al. is that the D/H 
in the local ISM is about 1.6 x 10"^ (Linsky et al. 1995, Piskunov et al. 
1997), while the relatively low metallicity of the Galaxy suggests that only 
a relatively modest fraction of the primordial D could have been destroyed 
(Tosi ct al. 1997). It thus seems that the lower D/H and correspondingly 
higher flh « O.lh^Q are more likely to be correct, although it is worrisome 
that the relatively high value Yp Ri 0.25 predicted by standard BBN for the 
primordial ^He abundance does not appear to be favored by the data (Olive 
et al. 1996, but cf. Sasselov &: Goldwirth 1995, Schramm & Turner 1997). 

White et al. (1993) have emphasized that X-ray observations of clusters, 
especially Coma, show that the abundance of baryons, mostly in the form 
of gas (which typically amounts to several times the mass of the cluster 
galaxies), is about 20% of the total cluster mass if h is as low as 0.5. For 
the Coma cluster they find that the baryon fraction within the Abell radius 
{1.5h-^ Mpc) is 

fb = > 0.009 + 0.050/^"^/^ (1.4) 

where the first term comes from the galaxies and the second from gas. If 
clusters are a fair sample of both baryons and dark matter, as they are 
expected to be based on simulations (Evrard, Metzler, &; Navarro 1996), 
then this is 2-3 times the amount of baryonic mass expected on the basis of 
BBN in an = 1, /t Ri 0.5 universe, though it is just what one would expect 
in a miivcrsc with ~ 0-3 (Steigman &; Felten 1995). The fair sample 
hypothesis implies that 

A recent review of X-ray measurements gas in a sample of clusters (White 
Sz Fabian 1995) finds that the baryon mass fraction within about 1 Mpc lies 
between 10 and 22% (for h = 0.5; the limits scale as /i~^/^), and argues 
that it is unlikely that (a) the gas could be clumped enough to lead to 
significant overestimates of the total gas mass — the main escape route 
considered in White et al. 1993 (cf. Gunn & Thomas 1996). The gas 
mass would also be overestimated if large tangled magnetic fields provide 
a significant part of the pressure in the central regions of some clusters 
(Loeb & Mao 1994, but cf. Felten 1996); this can be checked by observation 
of Faraday rotation of sources behind clusters (Kronberg 1994). If = 
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1, the alternatives are then either (b) that clusters have more mass than 
virial estimates based on the cluster galaxy velocities or estimates based 
on hydrostatic equilibrium (Balland & Blanchard 1995) of the gas at the 
measured X-ray temperature (which is surprising since they agree: Bahcall &; 
Lubin 1994), (c) that the usual BEN estimate of flj, is wrong, or (d) that the 
fair sample hypothesis is wrong (for which there is even some observational 
evidence: Loewenstein & Mushotzky 1996). It is interesting that there are 
indications from weak lensing that at least some clusters (e.g., for A2218 
see Squires et al. 1996; for this cluster the mass estimate from lensing 
becomes significantly higher than that from X-rays when the new ASCA 
satellite data, indicating that the temperature falls at large radii, is taken 
into account: Loewenstein 1996) may actually have extended halos of dark 
matter — something that is expected to a greater extent if the dark matter 
is a mixture of cold and hot components, since the hot component clusters 
less than the cold (Kofman et al. 1996). If so, the number density of clusters 
as a function of mass is higher than usually estimated, which has interesting 
cosmological implications (e.g. ag is higher than usually estimated). It 
is of course possible that the solution is some combination of alternatives 
(a)-(d). If none of the alternatives is right, then the only conclusion left is 
that no w 0.3. 

Notice that the rather high baryon fraction Jlj ~ 0.1(0.5//i)^ implied 
by the recent Tytler et al. measurements of low D/H helps resolve the 
cluster baryon crisis for J7 = 1 — it is escape route (c) above. With 
the higher fif, implied by the low D/H, there is now a "baryon cluster 
crisis" for Iow-JIq models! Even with a baryon fraction at the high end 
of observations, fh ^ 0.2(/i/0.5)~^/^, the fair sample hypothesis with this 
implies J^o ~ 0.5(/i/0.5)-V2. 

Another recent development is the measurement of the cluster baryon 
fraction using the Sunyaev-Zel'dovich effect (Myers et al. 1997), giving 

= (0.06ib0.01)/i~^. For h ~ 0.5, this is considerably lower than the X-ray 
estimates, and consistent with the Tytler et al. r^;, at the la level. The S-Z 
decrement is proportional to the electron density rig in the cluster, while the 
X-ray luminosity is proportional to n'^/. thus the S-Z measurement is likely to 
be less sensitive to small-scale clumping of the gas. If this first S-Z result for 
fh is confirmed by measurments on additional clusters, the cluster baryon 
data will have become an argument for $7 1. 



1.4-6 Cluster Morphology and Evolution 

Cluster Morphology. Richstone, Loeb, and Turner (1992) showed that 
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clusters are expected to be evolved — i.e. rather spherical and featureless 
— in low-il cosmologies, in which structures form at relatively high rcdshift, 
and that clusters should be more irregular in = 1 cosmologies, where they 
have formed relatively recently and are still undergoing significant merger 
activity. There are few known clusters that seem to be highly evolved 
and relaxed, and many that are irregular — some of which are obviously 
undergoing mergers now or have recently done so (see e.g. Burns et al. 1994). 
This disfavors low-fi models, but it remains to be seen just how low. Recent 
papers have addressed this. In one (Mohr et al. 1995) a total of 24 CDM 
simulations with = 1 or 0.2, the latter with f^A = or 0.8, were compared 
with data on a sample of 57 clusters. The conclusion was that clusters with 
the observed range of X-ray morphologies are very unlikely in the low-fi 
cosmologies. However, these simulations have been criticized because the 
J7o = 0.2 ones included rather a large amount of ordinary matter: $7^ = 0.1. 
(This is unrealistic both because h 0.8 provides the best fit for = 0.2 
CDM, but then the standard BBN upper limit is Qfe < 0.02/1-2 = 0.03; and 
also because observed clusters have a gas fraction of ~ 0.15(/i/0.5)-^/2.) 
Another study (Jing et al. 1995) using dissipationless simulations and not 
comparing directly to observational data found that ACDM with = 0.3 
and h = 0.75 produced clusters with some substructure, perhaps enough to 
be observationally acceptable (cf. Buote & Xu 1997). Clearly, this important 
issue deserves study with higher resolution hydrodynamic simulations, with 
a range of assumed Qi,, and possibly including at least some of the additional 
physics associated with the galaxies which must produce the metallicity 
observed in clusters, and perhaps some of the heat as well. Better statistics 
for comparing simulations to data may also be useful (Buote &; Tsai 1996). 
Cluster Evolution. There is evidence on the evolution of clusters at relatively 
low redshift, both in their X-ray properties (Henry et al. 1992, Castander 
et al. 1995, Ebeling et al. 1995) and in the properties of their galaxies. In 
particular, there is a strong increase in the fraction of blue galaxies with 
increasing redshift (the "Butcher-Oemler effect"), which may be difficult 
to explain in a low-density universe (Kauffmann 1994). Field galaxies do 
not appear to show such strong cvohition; indeed, a recent study concludes 
that over the redshift range 0.2 < z < 1.0 there is no significant evolution 
in the number density of "normal" galaxies (Steidel, Dickinson, &; Persson 
1994) . This is compatible with the predictions of various models, including 
CHDM with two neutrinos sharing a total mass of about 5 eV (see below), 
but the dependence of the number of clusters Ud on redshift can be a useful 
constraint on theories (Jing & Fang 1994, Bryan et al. 1994, Walter &: 
Klypin 1996, Eke et al. 1996). Some (e.g., Carlberg et al. 1997b; Bahcall, 
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Fan, & Cen 1997) have argued that presently available data show less fall off 
of nci{z) with increasing z than expected in J7 = 1 cosmologies, and already 
imply that < 1- These arguments are not yet entirely convincing because 
the cluster data at various redshifts are difficult to compare properly since 
they are rather inhomogeneous, and the data are not compared to a wide 
enough range of models (Gross et al. 1997). 



1.4-7 Early Structure Formation 

In linear theory, adiabatic density fluctuations grow linearly with the scale 
factor in an = 1 universe, but more slowly if 17 < 1 with or without a 
cosmological constant. As a result, if fluctuations of a certain size in an 

= 1 and an Qq = 0.3 theory are equal in amplitude at the present epoch 
{z = 0), then at higher redshift the fluctuations in the low-J7 model had 
higher amplitude. Thus, structures typically form earlier in low-O models 
than in = 1 models. 

Since quasars are seen at the highest redshifts, they have been used to 
try to constrain $7 = 1 theories, especially CHDM which because of the 
hot component has additional suppression of small-scale fluctuations that 
are presumably required to make early structure (e.g., Haehnelt 1993). The 
difficulty is that dissipationless simulations predict the number density of 
halos of a given mass as a function of redshift, but not enough is known 
about the nature of quasars — for example, the mass of the host galaxy 
— to allow a simple prediction of the number of quasars as a function of 
redshift in any given cosmological model. A more recent study (Katz et al. 
1994) concludes that very efficient cooling of the gas in early structures, and 
angular momentum transfer from it to the dark halo, allows for formation of 
at least the observed number of quasars even in models where most galaxy 
formation occurs late (cf. Eisenstein & Loeb 1995). 

Observers are now beginning to see significant numbers of what may be 
the central regions of galaxies in an early stage of their formation at redshifts 
z = 3 — 3.5 (Steidel et al. 1996; Giavalisco, Steidel, & Macchetto 1996) — 
although, as with quasars, a danger in using systems observed by emission 
is that they may not be typical. As additional observations (e.g., Lowenthal 
et al. 1996) clarify the nature of these objects, they can perhaps be used to 
constrain cosmological parameters and models. (This data is discussed in 
more detail in § 1.7.5| .) 



Another sort of high redshift object which may hold more promise for 
constraining theories is damped Lyman a systems (DLAS). DLAS are 
high column density clouds of neutral hydrogen, generally thought to be 
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protogalactic disks, which are observed as wide absorption features in quasar 
spectra (Wolfe 1993). They are relatively common, seen in roughly a third 
of all quasar spectra, so statistical inferences about DLAS are possible. At 
the highest redshift for which data was published in 1995, z = 3 — 3.4, 
the density of neutral gas in such systems in units of critical density was 
reported to be ^gas ~ 0.006, comparable to the total density of visible 
matter in the universe today (Lanzetta, Wolfe, & Turnshek 1995). Several 
papers (Mo & Miralda-Escude 1994, Kauffmann & Chariot 1994, Ma & 
Bertschinger 1994) pointed out that the CHDM model with fijy = 0.3 could 
not produce such a high ^gas- However, it has been shown that CHDM 
with Q.^ = 0.2 could do so (Klypin et al. 1995, cf. Ma 1995). The power 
spectrum on small scales is a very sensitive function of the total neutrino 
mass in CHDM models. This theory makes two crucial predictions: O^as 
must fall off at higher redshifts, and the DLAS at z ^ 3 mostly correspond 
to systems of internal rotation velocity or velocity dispersion less than about 
100 km s~^. This velocity can perhaps be inferred from the Doppler widths 
of the metal line systems associated with the DLAS. Preliminary reports 
regarding the amount of neutral hydrogen in such systems deduced from 
the latest data at redshifts above 3.5 appear to be consistent with the first 
of these predictions (Storrie-Lombardi et al. 1996). But a possible problem 
for the second (Wolfe 1996) is the large velocity widths and other statistical 
properties (Prochaska &: Wolfe 1997) of the metal line systems associated 
with the highest-redshift DLAS (e.g., Lu et al. 1996, at z = 4.4); if these 
actually indicate that a massive disk galaxy is already formed at such a high 
redshift, and if discovery of other such systems shows that they are not rare, 
that would certainly disfavor CHDM and other theories with relatively little 
power on small scales. However, other interpretations of such data which 
would not cause such problems for theories like CHDM are perhaps more 
plausible, since they are based on fairly realistic hydrodynamic simulations 
(Haehnelt, Steinmetz, & Ranch 1997; cf. § |1.7.5| ). More data will help 
resolve this question, along with DLAS models including dust absorption 
(Pei & Fall 1995), lensing (Bartelmann & Loeb 1996, Mailer et al. 1997), 
and effects of star formation (Kauffmann 1996, Somerville et al. 1997). 

One of the best ways of probing early structure formation would be to 
look at the main light output of the stars of the earliest galaxies, which is 
redshifted by the expansion of the universe to wavelengths beyond about 5 
microns today. Unfortunately, it is not possible to make such observations 
with existing telescopes; since the atmosphere blocks almost all such infrared 
radiation, what is required is a large infrared telescope in space. The Space 
Infrared Telescope Facility (SIRTF) has long been a high priority, and it 
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will be great to have access to the data such an instrument will produce 

when it is launched sometime in the next decade. In the meantime, the 
Near Infrared Canicra/Multi-Object Spectrograph (NICMOS), installed on 
Hubble Space Telescope in spring 1997, will help. Infrared spectrographs on 
the largest ground-baed telescopes will also be of great value. 

An alternative method is to look for the starlight from the earliest stars as 
extragalactic background infrared light (EBL). Although it is difficult to see 
this background light directly because our Galaxy is so bright in the near 
infrared, it may be possible to detect it indirectly through its absorption 
of TeV gamma rays (via the process 77 — e"^e~). Of the more than 
twenty active galactic nuclei (AGNs) that have been seen at ~ 10 GeV by 
the EGRET detector on the Compton Gamma Ray Observatory, only two 
of the nearest, Mk421 and Mk501, have also been clearly detected in TeV 
gamma rays by the Whipple Atmospheric Cerenkov Telescope (Quinn et 
al. 1996, Schubnell et al. 1996). Absorption of ~ TeV gamma rays from 
(AGNs) at redshifts z ~ 0.2 has been shown to be a sensitive probe of the 
EBL and thus of the era of galaxy formation (MacMinn &; Primack 1996; 
MacMinn, Somerville, &; Primack 1997). 



1.4 -8 Conclusions Regarding Q 

The main issue that has been addressed so far is the value of the cosmological 
density parameter J7. Arguments can be made for CIq 0.3 (and models 
such as ACDM) or for Q = 1 (for which the best class of models is probably 
CHDM), but it is too early to tell which is right. 

The evidence would favor a small « 0.3 if (1) the Hubble parameter 
actually has the high value Hq 75 favored by many observers, and the age 
of the universe to > 13 Gyr; or (2) the baryonic fraction = Mj^/Mtot in 
clusters is actually ~ 15%, about 3 times larger than expected for standard 
Big Bang Nucleosynthesis in an $7 = 1 universe. This assumes that standard 
BBN is actually right in predicting that the density of ordinary matter 
fife lies in the range 0.009 < Qbh"^ < 0.02. High-resolution high-redshift 
spectra are now providing important new data on primordial abundances 
of the light isotopes that should clarify the reliability of the BBN limits 
on rife- If the systematic errors in the ^He data are larger than currently 
estimated, then it may be wiser to use the deuterium upper limit < 
0.03, which is also consistent with the value n^/i^ 0.024 indicated by the 
only clear deuterium detection at high redshift, with the same D/H~ 2.4 x 
10~^ observed in two different low-metallicity quasar absorption systems 
(Tytler et al. 1996); this considerably lessens the discrepancy between 
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and fift. Another important constraint on fife will come from the new data 

on small angle CMB anisotropies — in particular, the location and height 
of the first Doppler (or acoustic, or Sakharov) peak (Dodclson, Gates, & 
Stebbins 1996; Jungman et al. 1996; Tegmark 1996), with the latest data 
consistent with low h 0.5—0.6 and high Qbh"^ ^ 0.025. The location of the 
first Doppler peak at angular wavenumber I ^ 250 indicated by the presently 
available data (Nettcrfield ct al. 1997, Scott ct al. 1996) is evidence in favor 
of a flat Tinivcrsc; ilo ^ 0.5 with A = is disfavored by this data (Lineweaver 
& Barbosa 1997). 

The evidence would favor O = 1 if (1) the POTENT analysis of galaxy 
peculiar velocity data is right, in particular regarding outflows from voids or 
the inability to obtain the present-epoch non-Gaussian density distribution 
from Gaussian initial fluctuations in a low-0 universe; or (2) the preliminary 
indication of high JIq and low Q,a from high-redshift Type la supernovae 
(Perlmutter et al. 1996) is confirmed. 

The statistics of gravitational lensing of quasars is incompatible with large 
cosmological constant A and low cosmological density 0,q. Discrimination 
between models may improve as additional examples of lensed quasars are 
searched for in large surveys such as the Sloan Digital Sky Survey. The 
era of structure formation is another important discriminant between these 
alternatives, low 0, favoring earlier structure formation, and = 1 favoring 
later formation with many clusters and larger-scale structures still forming 
today. A particularly critical test for models like CHDM is the evolution as 
a function of redshift of flgas in damped Lya systems. Reliable data on all 
of these issues is becoming available so rapidly today that there is reason to 
hope that a clear decision between these alternatives will be possible within 
the next few years. 

What if the data ends up supporting what appear to be contradictory 
possibilities, e.g. large JIq CLnd large Hq? Exotic initial conditions (e.g., 
"designer" primordial fluctuation spectra, cf. Hodges et al. 1990) or exotic 
dark matter particles beyond the simple "cold" vs. "hot" alternatives 
discussed in the next section (e.g., decaying 1-10 MeV tau neutrinos, 
Dodelson, Gyuk, & Turner 1994; volatile dark matter, Picrpaoli et al. 
1996) could increase the space of possible inflationary theories somewhat. 
But unless new observations, such as the new stellar parallaxes from the 
Hipparcos satellite, cause the estimates of Hq and to to be lowered, it 
may ultimately be necessary to go outside the framework of inflationary 
cosmological models and consider models with large scale spatial curvature, 
with a fairly large A as well as large f^o- This seems particularly unattractive, 
since in addition to implying that the universe is now entering a final 
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inflationary period, it means that inflation probably did not happen at 
the beginning of the universe, when it would solve the flatness, horizon, 
monopole, and structure-generation problems. Moreover, aside from the 
Hq — to problem, there is not a shred of reliable evidence in favor of A > 0, 
just increasingly stringent upper limits. Therefore, most cosmologists are 
rooting for the success of inflation-inspired cosmologies, with 0,q + = 1. 
With the new upper limits on A from gravitational lensing of quasars, 
number counts of elliptical galaxies, and high-redshift Type la supernovae, 
this means that the cosmological constant is probably too small to lengthen 
the age of the universe significantly. So one hopes that when the dust finally 
settles, Hq and to will both turn out to be low enough to be consistent 
with General Relativistic cosmology. But of course the universe is under no 
obligation to live up to our expectations. 



1.5 Dark Matter Particles 

1.5.1 Hot, Warm, and Cold Dark Matter 

The current limits on the total and baryonic cosmological density parameters 
have been summarized, and it was argued in particular that ~ 0.3 while 
Qb ^ 0.1. i^o > implies that the majority of the matter in the universe 
is not made of atoms. If the dark matter is not baryonic, what is it? 
Summarized here are the physical and astrophysical implications of three 
classes of elementary particle DM candidates, which are called hot, warm, 
and cold.! Table 1.3 gives a list of dark matter candidates, classified into 
these categories. 

Hot DM refers to low-mass neutral particles that were still in thermal 
equilibrium after the most recent phase transition in the hot early universe, 
the QCD confinement transition, which took place at Tqcd ~ 10^ MeV. 
Neutrinos are the standard example of hot dark matter, although other 
more exotic possibilities such as "majorons" have been discussed in the 
literature. Neutrinos have the virtue that z^p, v,,, and v^- are known to 



exist, and as summarized in § 1.5.3 there is experimental evidence that 



at least some of these neutrino species have mass, though the evidence is 
not yet really convincing. Hot DM particles have a cosmological number 
density roughly comparable to that of the microwave background photons, 
which implies an upper bound to their mass of a few tens of eV: m{u) = 
^uPo/^u = ^i/92h^ eV. Having iljy ~ 1 implies that free streaming destroys 

t Dick Bond suggested this terminology to me at the 1983 Moriond Conference, where I used 
it in my talk (Primack and Blumenthal, 1983). George Blumenthal and I had thought of this 
classification independently, but we used a more complicated terminology. 
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axion 

SUSY LSP ncutralino 
technibaryon 
pseudo Higgs 

shadow matter 
topological relics 
non-top. solitons 

Primordial BH 
jupiters 
brown dwarfs 
white dwarfs 
neutron stars 
stellar BH 
massive BH 



Weakly 
Interacting 
Massive 
Particles 



Massive 

Astrophysical 

Compact 

Halo 

Objects 



COLD 



' DARK 



MATTER 



neutrinos Ug (j^s?) I HOT 

majorons? | DARK MATTER 



gravitino 

right-handed v [ WARM DM 

decaying dark matter 

VOLATILE DM 



any adiabatic fluctuations smaller than supercluster size, ~ 10^^ M© (Bond, 
Efstathiou, & Silk 1980). With the COBE upper limit, HDM with adiabatic 
fluctuations would lead to hardly any structure formation at all, although 
Hot DM plus some sort of seeds, such as cosmic strings (see, e.g., Zanchin 
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et al. 1996), might still be viable. Another promising possibility is Cold + 
Hot DM with ~ 0.2 (CHDM, discussed in some detail below). 

Warm DM particles interact much more weakly than neutrinos. They 
decouple (i.e., their mean free path first exceeds the horizon size) at 
T ^ Tqcb, and they are not heated by the subsequent annihilation of 
hadronic species. Consequently their number density is expected to be 
roughly an order of magnitude lower, and their mass an order of magnitude 
higher, than hot DM particles. Fluctuations as small as large galaxy halos, 
> 10^^ Mq, could then survive free streaming. Pagels and Primack (1982) 
initially suggested that, in theories of local supersymmetry broken at ~ 10^ 
GeV, gravitinos could be DM of the warm variety. Other candidates have 
also been proposed, for example light right-handed neutrinos (Olive & 
Turner 1982). Warm dark matter does not lead to structure formation in 
agreement with observations, since the mass of the warm particle must be 
chosen rather small in order to have the power spectrum shape appropriate 
to fit observations such as the cluster autocorrelation function, but then it 
is too much like standard hot dark matter and there is far too little small 
scale structure (Colombi, Dodelson, & Widrow 1996). (This, and also the 
possibly promising combination of hot and more massive warm dark matter. 



will be discussed in more detail in S 1.7 below 



Cold DM consists of particles for which free streaming is of no cosmological 
importance. Two different sorts of cold DM consisting of elementary 
particles have been proposed, a cold Bose condensate such as axions, and 
heavy remnants of annihilation or decay such as super symmetric weakly 
interacting massive particles (WIMPs). As has been summarized above, a 
universe dominated by cold DM looks very much like the one astronomers 
actually observe, at least on galaxy to cluster scales. 



1.5.2 Cold Dark Matter Candidates 

The two sorts of particle candidates for cold dark matter that are best 
motivated remain supersymmetric Weakly Interacting Massive Particles 
(WIMPs) and the axion. They are both well motivated because both 
supersymmetry and Peccei-Quinn symmetry (associated with the axion) 
are key ideas in modern particle physics that were proposed independently 
of their implications for dark matter (for a review emphasizing direct and 
indirect methods of detecting both of these, see Primack, Seckel, &; Sadoulet 
1988). There are many other dark matter candidates whose motivations are 
more ad hoc (see Table 1.3) from the viewpoint of particle physics. But there 
is observational evidence that Massive Compact Halo Objects (MACHOs) 
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may comprise a substantial part of the mass of the Milky Way's dark matter 
halo. 

1.5.2.1 Axions 

Peccei-Quinn symmetry, with its associated particle the axion, remains 
the best solution known to the strong CP problem. A second-generation 

experiment is currently underway at LLNL (Hagmann et al. 1996) with 
sufficient sensitivity to have a chance of detecting the axions that might 
make up part of the dark matter in the halo of our galaxy, if the axion 
mass lies in the range 2-20 /xeV. However, it now appears that most of the 
axions would have been emitted from axionic strings (Battye & Shellard 
1994, 1997) and from the collapse of axionic domain walls (Nagasawa &: 
Kawasaki 1994), rather than arising as an axion condensate as envisioned 
in the original cosmological axion scenario. This implies that if the axion is 
the cold dark matter particle, the axion mass is probably ~ 0.1 meV, above 
the range of the LLNL experiment. While current experiments looking for 
either axion or supersymmetric WIMP cold dark matter have a chance of 
making discoveries, neither type is yet sufficiently sensitive to cover the full 
parameter space and thereby definitively rule out either theory if they do 
not detect anything. But in both cases this may be feasible in principle with 
more advanced experiments that may be possible in a few years. 

1.5.2.2 Supersymmetric WIMPs 

Prom the 1930s through the early 1970s, much of the development of 
quantum physics was a search for ever bigger symmetries, from spin and 
isospin to the Poincare group, and from electroweak symmetry to grand 
unified theories (GUTs). The larger the symmetry group, the wider the 
scope of the connections established between different elementary particles 
or other quantum states. The basic pattern of progress was to find the 
right Lie group and understand its role — SU(2) as the group connecting 
different states in the cases of spin and isospin; SU(3)xSU(2)xU(l) as 
the dynamical gauge symmetry group of the "Standard Model" of particle 
physics, connecting states without a gauge boson to states of the same 
particles including a gauge boson. Supersymmetry is a generalization of 
this idea of symmetry, since it mixes space-time symmetries, whose quantum 
numbers include the spin of elementary particles, with internal symmetries. 
It is based on a generalization of Lie algebra called graded Lie algebra, which 
involves anti-commutators as well as commutators of the operators that 
transform one particle state into another. Supersymmetry underlies almost 
all new ideas in particle physics since the mid-1970s, including superstrings. 
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Table 1.4. Supersymmetry 
A hypothetical symmetry between boson and fermion fields and interactions 



Spin 


Matter 


Forces 


Hypothetical 


Spin 




(fermions) 


(bosons) 


SuperpEirtners 




2 




graviton 


gravitino 


3/2 


1 




photon, W^,Z° 


photino, winos, zino, 


1/2 






gluons 


gluinos 




1/2 


quarks u,d,. . . 




squarks u,d,... 







leptons ejVe, ■ ■ ■ 




sleptons e, i'e, • • . 









Higgs bosons 


Higgsinos 


1/2 






axion 


axinos 





Note: Supersymmetric cold dark matter candidate particles are underlined. 



If valid, it is also bound to be relevant to cosmology. (Some reviews: Collins, 
Martin, & Squires 1989; de Boer 1994.) 

The simplest version of supersymmetry, which should be manifest at 
the GUT scale (~ 10^^ GeV) and below, has as its key prediction that 
for every kind of particle that we have learned about at the relatively 
low energies which even our largest particle accelerators can reach, there 
should be an as-yet-undiscovered "supersymmetric partner particle" with 
the same quantum numbers and interactions except that the spin of this 
hypothetical partner particle differs from that of the known particle by half 
a unit. For example, the partner of the photon (spin 1) is the "photino" 
(spin 1/2), and the partner of the electron (spin 1/2) is the "sclcctron" (spin 
0). Note that if a particle is a fermion (spin 1/2 or 3/2, obeying the Pauli 
exclusion principle), its partner particle is a boson (spin 0, 1,2). The familiar 
elementary particles of matter (quarks and leptons) arc all fermions, a fact 
that is responsible for the stability of matter, and the force particles are all 
bosons. Table 1.4 is a chart of the known families of elementary particles 
and their supersymmetric partners. It is these hypothetical partner particles 
among which we can search for the cold dark matter particle. The most 
interesting candidates are underlined. (As has already been mentioned, the 
gravitino is a warm dark matter particle candidate; this is discussed further 
below.) 

Note the parallel with Dirac's linking of special relativity and quantum 
mechanics in his equation for spin-1/2 particles (Griest 1996). In 
modern language, the resulting CPT invariance (under the combination of 
charge-conjugation C, replacing each particle with its antiparticle; parity 
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P, reversing the direction of each spatial coordinate; and time-reversal T) 
requires a doubling of the number of states: an anti-particle for every particle 
(except for particles, like the photon, which arc their own antiparticles). 

There are two other key features of supersymmetry that make it 
especially relevant to dark matter, i?-parity and the connection between 
supersymmetry breaking and the electroweak scale. The i?-parity of any 
particle is R = (^ — I^L+3B+2S ^ where L, B, and S are its lepton number, 
baryon number, and spin. Tims for an electron (L = 1, B = 0, S = 1/2) 
R = 1, and the same is true for a quark (L = 0, B = 1/3, S = 1/2) or a 
photon {L = 0, B = 0, S = 1). Indeed i? = 1 for all the known particles. 
But for a selectron {L = 1, B = 0, S = 1/2) or a photino {L = 0, B = 0, 
S = 1/2), the i?-parity is -1, or "odd". In most versions of supersymmetry, 
i?-parity is exactly conserved. This has the powerful consequence that the 
lightest i?-odd particle — often called the "lightest super symmetric partner" 
(LSP) — must be stable, for there is no lighter i?-odd particle for it to decay 
into. The LSP is thus a natural candidate to be the dark matter, as was first 
pointed out by Pagels &: me (1982), although as mentioned above the LSP 
in the early form of supersymmetry that we considered would have been a 
gravitino weighing about a keV, which would now be classified as warm dark 
matter. 

In the now-standard version of supersymmetry, there is an answer to the 
deep puzzle why there should be such a large difference in mass between the 
GUT scale Mqut ~ 10^^ GeV and the electroweak scale Mw = 80 GeV. 
Since both gauge symmetries are supposed to be broken by Higgs bosons 
which moreover must interact with each other, the natural expectation 
would be that Mqut ~ M\y. The super symmetric answer to this "gauge 
hierarchy" problem is that the masses of the weak bosons and all 
other light particles are zero until supersymmetry itself breaks. Thus, there 
is a close relationship between the masses of the supersymmetric partner 
particles and the electroweak scale. Since the abundance of the LSP is 
determined by its annihilation in the early universe, and the corresponding 
cross section involves exchanges of weak bosons or supersymmetric partner 
particles — all of which have electromagnetic-strength couplings and masses 
~ Mw — the cross sections will be o" ~ e^s/My^r (where s is the square 
of the center-of-mass energy) i.e., comparable to typical weak interactions. 
This in turn has the remarkable consequence that the resulting density 
of LSPs today corresponds to nearly critical density, i.e. f^LSP ~ 1- 
The LSP is typically a spin- 1/2 particle called a "neutralino" which is 
its own antiparticle — that is, it is a linear combination of the photino 
(supersymmetric partner of the photon), "zino" (partner of the weak 
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boson), "Higgsinos" (partners of the two Higgs bosons associated with 
electroweak symmetry breaking in supersymmetric theories), and "axinos" 
(partners of the axion, if it exists). In much of the parameter space, the 
neutrahno X is a "bino," a particular hnear combination of the photino and 
zino. Ah of these neutrahno LSPs are WIMPs, weakly interacting massive 
particles. Because of their large masses, several 10s to possibly 100s of GeV, 
these supersymmetric WIMPs would be dark matter of the "cold" variety. 

Having explained why supersymmetry is likely to be relevant to cold 
dark matter, one should also briefly summarize why supersymmetry is so 
popular with modern particle physicists. The reasons are that it is not only 
beautiful, it is even perhaps likely to be true. The supersymmetric pairing 
between bosons and fermions results in a cancellation of the high-energy (or 
"ultraviolet") divergences due to internal loops in Feynman diagrams. It 
is this cancellation that allows supersymmetry to solve the gauge hierarchy 
problem (how can Mqut/Mw be so big), and perhaps also unify gravity 
with the other forces ("superunification," "supergravity," "superstrings"). 
The one prediction of supersymmetry (Georgi, Quinn, &; Weinberg 1974) 
that has been verified so far is related to grand unification (Amaldi, de 
Boer, & Furstenau 1991). The way this is usually phrased today is that the 
three gauge couplings associated with the three parts of the standard model 
— the SU(3) "color" strong interactions, and the SU(2)xU(l) electroweak 
interactions — do not unify at any higher energy scale unless the effects of 
the supersymmetric partner particles are included in the calculation, and 
they do unify with the minimal set of partners (one partner for each of 
the known particles) as long as the partner particles all have masses not 
much higher than the electroweak scale M\y (which, as explained above, 
is expected if electroweak symmetry breaking is related to supersymmetry 
breaking). 

The expectations for the LSP neutrahno, including prospects for their 
detection in laboratory experiments and via cosmic rays, have recently been 
exhaustively reviewed (Jungman, Kamionkowski, &; Griest 1996). Several 
ambitious laboratory search experiments for LSPs in the mass range of tens 
to hundreds of GeV are now in progress (e.g., Shutt et al. 1996), and 
within the next few years they will have adequate sensitivity to probe a 
significant amount of the supersymmetric model parameter space. There 
are also hints of supersymmetric effects from recent experiments, which 
suggest that supersymmetry may be definitively detected in the near future 
as collider energy is increased — and also hint that the LSP may be rather 
light (Kane & Wells 1996), possibly even favoring the gravitino as the LSP 
(Dimopoulos et al. 1996). 
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1.5.2.3 MACHO s 

Meanwhile, the MACHO (Alcock et al. 1996a) and EROS (Ansari et 
al. 1996, Renault et al. 1996) experiments have detected microlensing of 
stars in the Large Magellanic Cloud (LMC). While the number of such 
microlensing events is small (six fairly convincing ones from two years of 
MACHO data discussed in their latest conference presentations, and one 
from three years of EROS observations), it is several times more than would 
be expected just from microlensing by the known stars. The MACHO data 
suggests that objects with a mass of O.S^qJM© are probably responsible 
for this microlensing, with their total density equal to ~ 20 — 50 percent 
of the mass of the Milky Way halo around ~ 20 kpc radius (Gates, Gyuk, 
Sz Turner 1996). Neither the EROS nor the MACHO groups have seen 
short duration microlensing events, which implies strong upper limits on 
the possible contribution to the halo of compact objects weighing less than 
about O.O5M0. While the MACHO masses are in the range expected for 
white dwarfs, there are strong observational limits (Flynn, Gould, & Bahcall 
1996) and theoretical arguments (Adams & Laughlin 1996) against white 
dwarfs being a significant fraction of the dark halo of our galaxy. Thus 
it remains mysterious what objects could be responsible for the observed 
microlensing toward the LMC. But the very large number of microlensing 
events observed toward the galactic bulge is probably explained by the 
presence of a bar aligned almost toward our position (Zhao, Rich, & Spergel 
1996; cf. Bissantz et al. 1996 for a dissenting view). Possibly the relatively 
small number of microlensing events toward the LMC represent lensing by 
a tidal tail of stars stretching toward us from the main body of the LMC 
(Zhao 1997); there is even some data on the colors and luminosities of stars 
toward the LMC suggesting that this may actually be true (D. Zaritsky, 
private communication 1997). 



1.5.3 Hot Dark Matter: Data on Neutrino Mass 

The upper limit on the electron neutrino mass is roughly 10-15 eV; the 

current Particle Data Book (Barnctt et al. 1996) notes that a more precise 
limit cannot be given since unexplained effects have resulted in significantly 
negative measurements of m(z>e)^ in recent precise tritium beta decay 
experiments. The (90% C.L.) upper limit on an effective Majorana neutrino 
mass 0.65 eV from the Heidelberg-Moscow ^^Ge neutrinoless double beta 
decay experiment (Balysh et al. 1995). The upper limits from accelerator 
experiments on the masses of the other neutrinos are m{u^) < 0.17 MeV 
(90% C.L.) and m(i/^) < 24 MeV (95% C.L.). Since stable neutrinos with 
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such large masses would certainly "overdose the universe" (i.e., prevent it 
from attaining its present age), the cosmological upper limits follow from the 
neutrino contribution to the cosmological density ilj^ = m,{v) / (^2h? cY) < 
Oq. There is a small window for an unstable with mass ~ 10 — 24 
MeV, which could have many astrophysical and cosmological consequences: 
relaxing the Big-Bang Nucleosynthesis bound on and N^,, allowing BBN 
to accommodate a low (less than 22%) primordial ^He mass fraction or 
high dcTiterium abundance, improving significantly the agreement between 
the CDM theory of structure formation and observations, and helping to 
explain how type II supernovae explode (Gyuk k, Turner 1995). 

But there is mounting astrophysical and laboratory data suggesting that 
neutrinos oscillate from one species to another, and therefore that they have 
non-zero mass. The implications if all these experimental results are taken 
at face value are summarized in Table 1.5. Of these experiments, the ones 
that are most relevant to neutrinos as hot dark matter are LSND and the 
higher energy Kamiokande atmospheric (cosmic ray) neutrinos. But the 
experimental results that are probably most secure are those concerning 
solar neutrinos, suggesting that some of the electron neutrinos undergo MSW 
oscillations to another species of neutrino as they travel through the sun (see, 
e.g., Hata & Langacker 1995, Bahcall 1996). 

The recent observation of events that appear to represent i>e 
oscillations followed hyr'e+P^n + e'^,n + p^D + j, with coincident 
detection of and the 2.2 MeV neutron-capture 7-ray in the Liquid 
Scintillator Neutrino Detector (LSND) experiment at Los Alamos suggests 
that Amg^ = |m(f^)^ - m(fe)^| > (Athanassopoulos et al. 1995, 1996). 
The analysis of the LSND data through 1995 strengthens the earlier LSND 
signal for — > oscillations. Comparison with exclusion plots from 
other experiments implies a lower limit Am^g = \m{i'^)'^ — m(fe)^| ^ 0.2 
eV^, implying in turn a lower limit nii, ^ 0.45 eV, or ^ 0.02(0.5//i)^. 
This implies that the contribution of hot dark matter to the cosmological 
density is larger than that of all the visible stars (Jl* f» 0.004 (Peebles 1993, 
eq. 5.150). More data and analysis are needed from LSND's ly^ Ve 
channel before the initial hint (Caldwell 1995) that Am^g « 6 eV^ can be 
confirmed. Fortunately the KARMEN experiment has just added shielding 
to decrease its background so that it can probe the same region of Am^g and 
mixing angle, with sensitivity as great as LSND's within about two years 
(Kleinfeller 1996). The Kamiokande data (Fukuda 1994) showing that the 
deficit of -E > 1.3 GeV atmospheric muon neutrinos increases with zenith 
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Table 1.5. Data Suggesting Neutrino Mass 

Solar Aml^ = 10"^ eV^, sin^26lox small 

Atm Ufj, deficit (9) Amf^y ~ lO'^ eV^ sin^ 26*^2, - 1 
Kamiokandc > 1.3 GeV 

Reactor i/g probably excludes y = e, so atm Vr or Vs 

BBN excludes v^^ — s- Vs with large mixing, so y = r 

LSND Am2^ « 1 - 10 eV^, sin226i^e small 
excludes x — fx, so solar —> Vs 

Cold + Hot Dark Matter Em^, w 5 eV /i|n 



angle suggests that t'^ ^ i^,- oscillations^ occur with an oscillation length 
comparable to the height of the atmosphere, implying that Am^^ ~ 10"'^ 
eV^ — which in turn implies that if either or u-j- have large enough 
mass ( ^ 1 eV) to be a hot dark matter particle, then they must be nearly 
degenerate in mass, i.e., the hot dark matter mass is shared between these 
two neutrino species. The much larger Super-Kamiokande detector is now 
operating, and we should know by about the end of 1996 whether the 
Kamiokande atmospheric neutrino data that suggested — > Vr oscillations 
will be confirmed and extended. Starting in 1997 there will be a long-baseline 
neutrino oscillation disappearance experiment to look for — > Vr with a 
beam of from the KEK accelerator directed at the Super-Kamiokande 
detector, with more powerful Fermilab-Soudan, KEK-Super-Kamiokande, 
and possibly CERN-Gran Sasso long-baseline experiments later. 

Evidence for non-zero neutrino mass evidently favors CHDM, but it 
also disfavors low-Q models. Because free streaming of the neutrinos 
damps small-scale fluctuations, even a little hot dark matter causes reduced 
fluctuation power on small scales and requires substantial cold dark matter 
to compensate; thus evidence for even 2 eV of neutrino mass favors large ft 
and would be incompatible with a cold dark matter density Jlc as small as 0.3 

t The Kamiokande data is consistent with atmospheric Ufj, oscillating to any other neutrino species 
y with a large mixing angle O/^y. But as summarized in Table 1.5 (see further discussion and 
references in, e.g., Primack et al. 1995, hereafter PHKC95; Fuller, Primack, & Qian 1995) U/j, 
oscillating to We with a large mixing angle is probably inconsistent with reactor and other data, 
and Ufj, oscillating to a sterile neutrino Us (i.e., one that does not interact via the usual weak 
interactions) with a large mixing angle is inconsistent with the usual Big Bang Nucleosynthesis 
constraints. Thus, by a process of elimination, if the Kamiokande data indicating atmospheric 
neutrino oscillations is right, the oscillation is Ufj, —* u-r- 
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(PHKC95). Allowing and the tilt to vary, CHDM can fit observations 
over a somewhat wider range of values of the Hubble parameter h than 
standard or tilted CDM (Pogosyan & Starobinsky 1995a, Liddle et al. 
1996b). This is especially true if the neutrino mass is shared between 
two or three neutrino species (Holtzman 1989; Holtzman Sz Primack 1993; 
PHKC95; Pogosyan & Starobinsky 1995b; Babu, Schaefer, & Shafi 1996), 
since then the lower neutrino mass results in a larger free-streaming scale 
over which the power is lowered compared to CDM. The result is that 
the cluster abundance predicted with Oj, « 0.2 and h « 0.5 and COBE 
normalization (corresponding to erg 0.7) is in reasonable agreement with 
observations without the need to tilt the model (Borgani et al. 1996) and 
thereby reduce the small-scale power further. (In CHDM with a given fl^, 
shared between Ni, = 2 or 3 neutrino species, the linear power spectra are 
identical on large and small scales to the Ni, = 1 case; the only difference is 
on the cluster scale, where the power is reduced by ~ 20% (Holtzman 1989, 
PHKC95, Pogosyan & Starobinsky 1995). 

1.6 Origin of Fluctuations: Inflation and Topological Defects 

1.6.1 Topological defects 

A fundamental scalar field, the Higgs field, is invoked by particle theorists 
to account for the generation of mass; one of the main goals of the next 
generation of particle accelerators, including the Large Hadron Collider at 
CERN, will be to verify the Higgs theory for the generation of the mass of 
the weak vector bosons and all the lighter elementary particles. Another 
scalar field is required to produce the vacuum energy which may drive 
cosmic inflation (discussed in the next section). Scalar fields can also create 
topological defects that might be of great importance in cosmology. The 
basic idea is that some symmetry is broken wherever a given scalar fleld 
(j) has a non-vanishing value, so the dimensionality of the corresponding 
topological defect depends on the number of components of the scalar field: 
for a single-component real scalar field, ^(r) = defines a two-dimensional 
surface in three-dimensional space, a domain wall; for a complex scalar field, 
the real and imaginary parts of 0(r) = define a one-dimensional locus, a 
cosmic string; for a three-component (e.g., isovector) field, 4>i{r) = for 
i = 1,2,3 is satisfied at isolated points, monopoles; for more than three 
components, one gets textures that are not topologically stable but which 
can seed structure in the universe as they unwind. 

To see how this works in more detail, consider a cosmic string. For 
the underlying field theory to permit cosmic strings, we need to couple 
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a complex scalar field to a single-component (i.e., U(l)) gauge field 
Aa, like the electromagnetic field, in the usual way via the substitution 
da Da = {da — icAa), SO that the scalar field derivative term in the 
Lagrangian becomes /Z/j^ = {Datpl"^- Then if the scalar field (p gets a non-zero 
value by the usual Higgs "spontaneous symmetry breaking" mechanism, the 
gauge symmetry is broken because the field has a definite complex phase. 
But along a string where (f) = the symmetry is restored. As one circles 
around the string at any point on it, the complex phase of (j){f) in general 
makes one, or possibly n > 1, complete circles 2mT. But since such a 
phase rotation can be removed at large distance from the string by a gauge 
transformation of cj) and , the energy density associated with this behavior 
of (j) C,D4> ^ at large distances, and therefore the energy ^ per unit length 
of string is finite. Since it would require an infinite amount of energy to 
unwind the phase of 4> at infinity, however, the string is topologically stable. 
If the field theory describing the early universe includes a U(l) gauge field 
and associated complex Higgs field (j), a rather high density of such cosmic 
strings will form when the string field (p acquires its nonzero value and breaks 
the U(l) symmetry. This happens because there is no way for the phase of cf) 
to be aligned in causally disconnected regions, and it is geometrically fairly 
likely that the phases will actually wrap around as required for a string to go 
through a given region (Kibble 1976). The string network will then evolve 
and can help cause formation of structure after the universe becomes matter 
dominated, as long as the string density is not diluted by a subsequent period 
of cosmic inflation (on the difficult problem of combining cosmic defects 
and inflation, see, e.g., Hodges & Primack 1991). A similar discussion can 
be given for domain walls and local (gauge) monopoles, but these objects 
are cosmologically pathological since they dominate the energy density and 
"overdose" the universe. But cosmic strings, a sufficiently low density of 
global (i.e., non-gauged monopoles), and global textures are potentially 
interesting for cosmology (recent reviews include Vilenkin & Shellard 1994, 
Hindmarsh & Kibble 1995, Shellard 1996). Cosmic defects are the most 
important class of models producing non-Gaussian fluctuations which could 
seed cosmic structure formation. Since they arc geometrically extended 
objects, they correspond to non-local non-Gaussian fluctuations (Kofman et 
al. 1991). 

The parameter usually quoted in the dimensionless form G/Lt (where G is 
Newton's constant), is the key parameter of the theory of cosmic strings. The 
value required for the COBE normalization is Gixq = GfixlO^ = 1 — 2 (recent 
determinations include GfiQ = 1.7ib0.7, Perivolaropoulos 1994; 2, Coulson et 
al. 1994; {l.ObtoH), ^Hen et al. 1996; 1.7, Allen et al. 1997). This is close 
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enough to the value required for structure formation, G/j, = (2.2 — 2.8)6g x 
10~^ (Albrecht & Stebbins 1992), with the smaller value for cosmic strings 
plus cold dark matter and the higher value for cosmic strings plus hot dark 
matter, so that the necessary value of the biasing factor bs is 1.3-3, which 
is high (probably leading to underproduction of clusters, and large-scale 
velocities that are low compared to observations — cf. Perivolaropoulos & 
Vachaspati 1994), but perhaps not completely crazy. (Here feg is the factor 
by which galaxies must be more clustered than dark matter, on a scale 
of 8h^^ Mpc.) Since generically G/i ~ {M/mpi)'^, where M is the energy 
scale at which the string field ^ acquires its nonzero value, the fact that 
G/j, ~ 10~^, corresponding to M at roughly the Grand Unification scale, 
is usually regarded as a plus for the cosmic string scenario. (Even though 
there is no particular necessity for cosmic strings in GUT scenarios, GUT 
groups larger than the minimal SU(5) typically do contain the needed extra 
U(l)s.) Moreover, the required normalization is well below the upper limit 
obtained from the requirement that the gravitational radiation generated by 
the evolution of the string network not disrupt Big Bang Nucleosynthesis, 
Gfi ^ 6 X 10~^. However, there is currently controversy whether it is also 
below the upper limit from pulsar timing, which has been determined to be 
G// < 6 X 10-'^ (Thorsett & Dewey 1996) vs. G// < 5 x 10"^ (McHugh et al. 
1996; cf. Caldwell, Battye, & Shellard 1996). 

As for cosmic strings, the COBE normalization for global texture models 
also implies a high bias bs ~ 3.4 for h = 0.7 (Bennett & Rhie 1993), although 
the needed bias is somewhat lower for Q ~ 0.3 (Pen & Spergel 1995). The 
latest global defect simulations (Pen, Seljak, &: Turok, 1997) show that the 
matter power spectrum in all such models also has a shape very different 
than that suggested by the available data on galaxies and clusters. 

But both cosmic string and global defect models have a problem which 
may be even more serious: they predict a small- angle CMB fluctuation 
spectrum in which the first peak is at rather high angular wavenumber 
i ^ 400 (Crittenden & Turok 1995, Durrer et al. 1996, Magueijo et 
al. 1996) and in any case is rather low in amplitude, partly because 
of incoherent addition of scalar, vector, and tensor modes, according to 
the latest simulations (strings: Allen et al. 1997; global defects: Pen, 
Seljak, k Turok 1997; cf. Albrecht, Battye, & Robinson 1997). This is 
in conflict with the currently available small-angle CMB data (Netterfield 
et al. 1997, Scott et al. 1996), which shows a peak at ^ ~ 250 and a 
drop at £ ^ 400, as predicted by flat (f^o -|- = 1) CDM-type models. 
Since the small-angle CMB data is still rather preliminary, it is premature 
to regard the cosmic defect models as being definitively ruled out. It 



Dark Matter and Structure Format- on 



55 



will be interesting to see the nature of the predicted galaxy distribution 
and CMB anisotropies when more complete simulations of cosmic defect 
models are run. This is more difficult than simulating models with the 
usual inflationary fluctuations, both because it is necessary to evolve the 
defects, and also because the fact that these defects represent rare but high 
amplitude fluctuations necessitates a careful treatment of their local effects 
on the ordinary and dark matter. It may be difficult to sustain the effort such 
calculations require, because the poor agreement between the latest defect 
simulations and current small-angle CMB data does not bode well for defect 
theories. Fortunately, there have been significant technical breakthroughs 
in calculational techniques (cf. Allen et al. 1997, Pen et al. 1997). 

1.6.2 Cosmic Inflation: Introduction 

The basic idea of inflation is that before the universe entered the present 
adiabatically expanding Friedmann era, it underwent a period of de Sitter 
exponential expansion of the scale factor, termed inflation (Guth 1981). 
Actually, inflation is never precisely de Sitter, and any superluminal 
(faster-than-light) expansion is now called inflation. Inflation was originally 
invented to solve the problem of too many GUT monopoles, which, as 
mentioned in the previous section, would otherwise be disastrous for 
cosmology. 

The de Sitter cosmology corresponds to the solution of Friedmann's 
equation in an empty universe (i.e., with p = 0) with vanishing curvature 
(k = 0) and positive cosmological constant (A > 0). The solution is 
a = ttoC^^, with constant Hubble parameter H = (A/3)^/^. There are 
analogous solutions for A; = +1 and k = —1 with a oc cosh Ht and 
a oc sinh Ht respectively. The scale factor expands exponentially because the 
positive cosmological constant corresponds effectively to a negative pressure, 
de Sitter space is discussed in textbooks on general relativity (for example, 
Rindler 1977, Hawking & Ellis 1973) mainly for its geometrical interest. 
Until cosmological inflation was considered, the chief significance of the de 
Sitter solution in cosmology was that it is a limit to which all indefinitely 
expanding models with A > must tend, since as a ^ oo, the cosmological 
constant term ultimately dominates the right hand side of the Friedmann 
equation. 

As Guth (1981) emphasized, the de Sitter solution might also have been 
important in the very early universe because the vacuum energy that plays 
such an important role in spontaneously broken gauge theories also acts as 
an effective cosmological constant. A period of de Sitter inflation preceding 
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Table 1.6. Inflation Summary 



Problem Solved 

Horizon Homogeneity, Isotropy, Uniform T 

Flatness/ Age Expansion and gravity balance 

"Dragons" Monopoles, doman walls,. . . banished 

Structure Small fluctuations to evolve into galaxies, 

clusters, voids 



Cosmological constant A > space repels space, so the more space the more 
repulsion, de Sitter exponential expansion a oc e^*. 

Inflation is exponentially accelerating expansion caused by efl^ective cosmological 
constant ( "false vacuum" energy) associated with hypothetical scalar field 
("inflaton"). 



Known 

Goal of LHC 
Early universe 



Forces of Nature Spin 

Gravity 2 

Strong, weak, and electromagnetic 1 

Mass (Higgs Boson) 

Inflation (Inflaton) 



Inflation lasting only ^10~ s sufiices to solve all the problems listed above. 
Universe must then convert to ordinary expansion through conversion of false to 
true vacuum ("re-" heating). 



ordinary radiation-dominated Priedmann expansion could explain several 

features of the observed universe that otherwise appear to require very 
special initial conditions: the horizon, flatness/age, monopole, and structure 
formation problems. (See Table 1.6.) 

Let us illustrate how inflation can help with the horizon problem. At 
recombination (p+ + ^ H), which occurs at a/ag ~ 10~^, the 
mass encompassed by the horizon was Mh ~ IO^^Mq, compared to 
Mh,o ~ 10^^ M0 today. Equivalently, the angular size today of the causally 
connected regions at recombination is only AO ~ 3°. Yet the fluctuation 
in temperature of the cosmic background radiation from different regions 
is very small: AT/T ~ 10~^. How could regions far out of causal 
contact have come to temperatures that are so precisely eqiial? This is 
the "horizon problem". With inflation, it is no problem because the entire 
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observable universe initially lay inside a single causally connected region that 
subsequently inflated to a gigantic scale. Similarly, inflation exponentially 
dilutes any proceeding density of monopoles or other unwanted relics (a 
modern version of the "dragons" that decorated the unexplored borders of 
old maps). 

In the first inflationary models, the dynamics of the very early universe 

was typically controlled by the self-energy of the Higgs field associated with 
the breaking of a Grand Unified Theory (GUT) into the standard 3-2-1 
model: GUT^ SU {3)coior ® [SU (2) U {l)]eiectroweak- This occurs when the 
cosmological temperature drops to the unification scale Tqut ~ 10^^ GeV at 
about 10~^^ s after the Big Bang. Guth (1981) initially considered a scheme 
in which inflation occurs while the universe is trapped in an unstable state 
(with the GUT unbroken) on the wrong side of a maximum in the Higgs 
potential. This turns out not to work: the transition from a de Sitter to a 
Priedmann universe never finishes (Guth Sz Weinberg 1981). The solution in 
the "new inflation" scheme (Linde 1982; Albrecht and Steinhardt 1982) is for 
inflation to occur after barrier penetration (if any). It is necessary that the 
potential of the scalar field controlling inflation ("in/Zaton") be nearly flat 
(i.e., decrease very slowly with increasing inflaton field) for the infiationary 
period to last long enough. This nearly fiat part of the potential must then 
be followed by a very steep minimum, in order that the energy contained in 
the Higgs potential be rapidly shared with the other degrees of freedom 
("reheating"). A more general approach, "chaotic" infiation, has been 
worked out by Linde (1983, 1990) and others; this works for a wide range of 
inflationary potentials, including simple power laws such as A0^. However, 
for the amplitude of the fluctuations to be small enough for consistency with 
observations, it is necessary that the inflaton self-coupling be very small, for 
example A ~ 10~^^ for the 0^ model. This requirement prevents a Higgs field 
from being the inflaton, since Higgs fields by definition have gauge couplings 
to the gauge field (which are expected to be of order unity), and these would 
generate self-couplings of similar magnitude even if none were present. Both 
the Higgs and inflaton are hypothetical fundamental (or possibly composite) 
scalar flelds (see Table 1.6). 

It turns out to be necessary to inflate by a factor ^ e^^ in order to solve 
the flatness problem, i.e., that ~ 1- (With ~ 10~^^ s during the de 
Sitter phase, this implies that the inflationary period needs to last for only 
a relatively small time r ^ 10~^^ s.) The "flatness problem" is essentially 
the question why the universe did not become curvature dominated long 
ago. Neglecting the cosmological constant on the assumption that it is 
unimportant after the inflationary epoch, the Friedmann equation can be 
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written 

a) = ^30^(^)^ - W ^ 

where the first term on the right hand side is the contribution of the energy 
density in relativistic particles and g{T) is the effective number of degrees of 
freedom. The second term on the right hand side is the curvature term. Since 
aT ~ constant for adiabatic expansion, it is clear that as the temperature 
T drops, the curvature term becomes increasingly important. The quantity 
K = k/{aT)'^ is a dimensionless measure of the curvature. Today, \K\ = 
\i^ — l\H^/T^ < 2 X 10~^^. Unless the curvature exactly vanishes, the 
most "natural" value for K is perhaps K 1. Since inflation increases 
a by a tremendous factor e^'^ at essentially constant T (after reheating), 
it increases aT by the same tremendous factor and thereby decreases the 
curvature by that factor squared. Setting e"^^"^ 2 x 10~^^ gives the needed 
amount of inflation: Ht ^ 66. This much inflation turns out to be enough 
to take care of the other cosmological problems mentioned above as well. 

Of course, this is only the minimum amount of inflation needed; the actual 
inflation might have been much greater. Indeed it is frequently argued that 
since the amount of inflation is a tremendously sensitive function of (e.g.) 
the initial value of the inflaton fleld, it is extremely likely that there was 
much more inflation than the minimum necessary to account for the fact 
that the universe is observed to be nearly flat today. It then follows (in the 
absence of a cosmological constant today) that = 1 to a very high degree 
of accuracy. A way of evading this that has recently been worked out is 



discussed in § 1.6.6 



1.6.3 Inflation and the Origin of Fluctuations 

Thus far, it has been sketched how inflation stretches, flattens, and smooths 
out the universe, thus greatly increasing the domain of initial conditions 
that could correspond to the universe that we observe today. But inflation 
also can explain the origin of the fluctuations necessary in the gravitional 
instability picture of galaxy and cluster formation. Recall that the very 
existence of these fluctuations is a problem in the standard Big Bang picture, 
since these fluctuations are much larger than the horizon at early times. How 
could they have arisen? 

The answer in the inflationary universe scenario is that they arise from 
quantum fluctuations in the inflaton fleld (f) whose vacuum energy drives 
inflation. The scalar fluctuations 5(f) during the de Sitter phase are of the 
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order of the Hawking temperature H/2TT. Because of these fluctuations, 
there is a time spread At ~ S(j)/(f) during which different regions of the 
same size complete the transition to the Fricdmann phase. The result is 
that the density fluctuations when a region of a particular size re-enters 
the horizon are equal to (Guth & Pi 1982; see Linde 1990 for alternative 
approaches) 6h = {Sp/p)h ~ At/iff = HAt. The time spread At can be 
estimated from the equation of motion of (the free Klein-Gordon equation 
in an expanding universe): + 3//<^ = —(dV/dcf)). Neglecting the (f) term, 
since the scalar potential V must be very flat in order for enough inflation 
to occur (this is called the "slow roll" approximation), Ri —V'/{3H), so 
5h ~ H^/V ~ y^/^/y'. Unless there is a special feature in the potential 
y((/)) as (p rolls through the scales of importance in cosmology (producing 
such "designer inflation" features generally requires fine tuning — see e.g. 
Hodges et al. 1990), V and V will hardly vary there and hence 6h will 
be essentially constant. These are fluctuations of all the contents of the 
universe, so they are adiabatic fluctuations. 

Thus inflationary models typically predict a nearly constant curvature 
spectrum 6h = constant of adiabatic fluctuations. Some time ago Harrison 
(1970), Zel'dovich (1972), and others had emphasized that this is the only 
scale-invariant (i.e., power-law) fluctuation spectrum that avoids trouble at 
both large and small scales. If 6h oc M^", where Mjj is the mass inside 
the horizon, then if —a is too large the universe will be less homogeneous 
on large than small scales, contrary to observation; and if a is too large, 
fluctuations on sufficiently small scales will enter the horizon with Su ^ ^ 
and collapse to black holes (see e.g. Carr, Gilbert, & Lidsey 1995, Bullock 
& Primack 1996); thus a « 0. The a = case has come to be known as the 
Zel'dovich spectrum. 

Inflation predicts more: it allows the calculation of the value of the 
constant 5h in terms of the properties of the scalar potential V^cj)). 
Indeed, this proved to be embarrassing, at least initially, since the 
Coleman- Weinberg potential, the first potential studied in the context of 
the new inflation scenario, results in Sh ~ 10^ (Guth Sz Pi 1982) some 
six orders of magnitude too large. But this does not seem to be an 
insurmountable difficulty; as was mentioned above, chaotic inflation works, 
with a sufficiently small self-coupling. Thus inflation at present appears to 
be a plausible solution to the problem of providing reasonable cosmological 
initial conditions (although it sheds no light at all on the fundamental 
question why the cosmological constant is so small now). Many variations of 
the basic idea of inflation have been worked out, and the following sections 
will discuss two recent developments in a little more detail. Linde (1995) 
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Table 1.7. Linde's Classification of Inflation Models 



How Inflation Begins 

Old Inflation Tinitiai high, 4>in « is false vacuum until phase transition 
Ends by bubble creation; Reheat by bubble collisions 

New Inflation Slow roll down V{(l)), no phase transition 

Chaotic Inflation Similar to New Inflation, but essentially arbitrary: 
any region with ^(p"^ + ^{diCp)'^ ^ ^(0) inflates 

Extended Inflation Like Old Inflation, but slower (e.g., power aoct^), 
so phase transition can finish 

Potential y(^) During Inflation 

Chaotic typically V{^) = can also use V = Voe°"^, etc. 

a oc iP, p= IGw/a'^ > 1 

How Inflation Ends 

First-order phase transition — e.g.. Old or Extended inflation 
Faster rolling — > oscillation — e.g.. Chaotic y((/))^A(/)" 
"Waterfall" — rapid roll of a triggered by slow roll of (f) 

(Re) heating 

Decay of inflatons 

"Preheating" by parametric resonance, then decay 

Before Inflation? 

Eternal Inflation? Can be caused by 

• Quantum Stf) ~ H/2Tr > rolhng = (j)At = (pH''^ « V'/V 

• Mouopol(-s or ()tli(>r tojiolop^ical defects 



recently classified these inflationary models in an interesting and useful way: 
see Table 1.7. 



1.6.4 Eternal Inflation 

Vilenkin (1983) and Linde (1986, 1990) pointed out that if one extrapolates 
inflation backward to try to imagine what might have preceeded it, in many 
versions of inflation the answer is "eternal inflation" : in most of the volume 
of the universe inflation is still happening, and our part of the expanding 
universe (a region encompassing far more than our entire cosmic horizon) 
arose from a tiny part of such a region. To see how eternal inflation works, 



Dark Matter and Structure Format-on 



61 



consider the simple chaotic model with V{(j)) = [m? /2)4p'. During the de 
Sitter Hubble time , where as usual H'^ = {SttG /'i)V , the slow rolling 
of (j) down the potential will reduce it by 

A^^m=--At = ^. (1.7) 

Here mpi is the Planck mass (see Table 1.1). But there will also be quantum 
fluctuations that will change (j) up or down by 

dcP=^= (1.8) 

These will be equal for = rri^Ji^ /2m}/'^ , y{<P*) = {m/^mpi)m'^pi. If(/) ^ 0*, 
positive quantum fluctuations dominate the evolution: after At ~ , an 
initial region becomes ~ regions of size ~ , in half of which (j) increases 
to (f) + 6(j). Since H oc <j), this drives inflation faster in these regions. Various 
mechanisms probably cut this off as (/) — > m?pi/m and V — > rrip^ — for further 
discussion and references, see Linde (1995). Thus, although (f) at any given 
point is likely eventually to roll down the potential and end inflation, in 
most of the volume of the metauniverse <j) > c/)^ and inflation is proceeding 
at a very fast rate. 

Eternal Inflation is eternal in the sense that, once started, it never ends. 
But it remains uncertain whether or not it could have begun an infinite 
length of time ago. Assuming the "weak energy condition" T^i,V^V'^ > for 
all timelike vectors V^, i.e. that any observer will measure a positive energy 
density, Borde & Vilenkin (1994) proved that a future-eternal inflationary 
model cannot be extended into the infinite past. However, Borde &: Vilenkin 
(1997) have recently shown that the weak energy condition is quite likely to 
be violated in inflating spacetimes (except the open universe inflation models 
discussed below, § 1.6.6| ), so a "steady-state" eternally inflating universe may 
be possible after all, with no beginning as well as no end. 



1.6.5 A Supersymmetric Inflation Model 

We have already considered, in connection with cold dark matter candidates, 
why supersymmetry is likely to be a feature of the fundamental theory of the 
particle interactions, of which the present "Standard Model" is presumably 
just a low-energy approximation. If the higher-energy regime within which 
cosmological inflation occurs is described by a supersymmetric theory, there 
are new cosmological problems that initially seemed insuperable. But recent 
work has suggested that these problems can plausibly be overcome, and 
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that supersymmetric inflation might also avoid the fine-tuning otherwise 
required to explain the small inflaton coupling corresponding to the COBE 
fluctuation amplitude. Here the problems will be briefly summarized, 
and an explanation will be given of how one such model, due to Ross &; 
Sarkar (1996; hereafter RS96) overcomes them. (An interesting alternative 
supersymmetric approach to inflation is sketched in Dine et al. 1996.) 

When Pagels and I (1982) flrst suggested that the lightest supersymmetric 
partner particle (LSP), stable because of R-parity, might be the dark 
matter particle, that particle was the gravitino in the early version of 
supersymmetry then in fashion. Weinberg (1982) immediately pointed 
out that if the gravitino were not the LSP, it could be a source or real 
trouble because of its long lifetime ~ Mp^/m^^^ ~ ("1-3/2 /TeV)~^ 10^ 
s, a consequence of its gravitational-strength coupling to other fields. 
Subsequently, it was realized that supersymmetric theories can naturally 
solve the gauge hierarchy problem, explaining why the electroweak scale 
Mew ~ 10^ GeV is so much smaller than the GUT or Planck scales. In 
this version of supersymmetry, which has now become the standard one, 
the gravitino mass will typically be 7713/2 ~ TeV; and the late decay of 
even a relatively small number of such massive particles can wreck BBN 
and/or the thermal spectrum of the CBR. The only way to prevent this is 
to make sure that the reheating temperature after inflation is sufficiently 
low: Trh ^ 2 X 10^ GeV (for 7773/2 = TeV) (Ellis, Kim, & Nanopoulos 1984, 
Ellis et al. 1992). 

This can be realized in supergravity theories rather naturally (RS96). 
Define M = Mpi/(87r)V2 = 2.4 x 10^^ Gev. Break GUT by the Higgs 
field X with vacuum expectation value (vev) < x >~ 10^^ GeV. Break 
supersymmetry by a gaugino condensate < AA >~ (lO^^GeV)^; then the 
gravitino mass is 7773/2 ~< AA > /M^ ~ TeV. Inflation is expected to inhibit 
such breaking, so it must occur afterward. The inflaton superpotential has 
the form / = A^M/((/>/M), with the corresponding potential 



with minimum at 00 . Demanding that at this minimum the potential 
actually vanishes V{<po) = 0, i.e., that the cosmological constant vanishes, 
implies that /((/>o) = {dl /d(p),f,f^ = 0. The simplest possibility is J = 
A'^ ((j) — (j)Q)'^ / M . Requiring that dV/dcplo = for a sufficiently fiat potential, 
implies that (pQ = M and that the second derivative also vanishes at the 




(1.9) 
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origin; thus 




(1.10) 



(Holman, Raymond, &: Ross 1984). This particular inflaton potential is of 
the "new inflation" type, and corresponds to tilt Up = 0.92 and a number of 
e-folds during inflation 



assuming that the starting value of the inflaton field (pin is suflficiently close 
to the origin (which has relatively small but nonvanishing probability — the 
(j) field presumably has a broad initial distribution). Matching The COBE 
fluctuation amplitude requires that A/M = 1.4 x 10"'^, which in turn implies 
that ~ 10^, ~ A^/M ~ 10^^ GeV, Trh ~ 10^ GeV (parametric 
resonance reheating does not occur). Such a low reheat temperature insures 
that there will be no gravitino problem, and requires that the baryon 
asymmetry be generated by clcctrowcak baryogenesis — which appears to 
be viable as long as the theory contains adequately large CP violation. 

Note the following features of the above scenario: inflation occurs at an 
energy scale far below the GUT scale, so there is essentially no gravity 
wave contribution to the large-angle CMB fluctuations (i.e., T/S ^ 0) even 
through there is significant tilt (rip = 0.92 for the particular potential above); 
there is a low reheat temperature, so electroweak baryogenesis is required; 
and the universe is predicted to be very flat since there are many more 
e-folds than required to solve the flatness problem. 



Can inflation produce a region of negative curvature larger than our present 
horizon — for example, a region with f^o < 1 and A = 0? The old approach 
to this problem was to imagine that there might be just enough inflation 
to solve the horizon problem, but not quite enough to over solve the flatness 
problem, e.g. A'^ ~ 60 (Steinhardt 1990). This requires fine tuning, but 
the real problem with this approach is that the resulting region will not be 
smooth enough to agree with the small size of the quadrupole anisotropy Q 
measured by COBE. According to the Grischuk-Zel'dovich (1978) theorem 
(cf. Garcia-Bellido et al. 1995), 5 ~ 1 fluctuations on a super- horizon scale 
L > Hq^ imply Q ~ (Li^o)"^- COBE measured Qrms < 2 x 10"^, which 
implies in turn that the region containing our horizon must be homogeneous 
on a scale L > 500Hq^, i.e. N > 70, |1 - f^ol ^ lO"! 




(1.11) 



1.6.6 Inflation with < 1 



64 



Dark Matter and Structure Formation 



Recently a new approach was discovered, based on the fact that a bubble 
created from de Sitter space by quantum tunneling tends to be spherical and 
homogeneous if the tunneling is sufficiently improbable. The interior of such 
bubbles are quite empty, i.e., they are a region of negative curvature with 

^ 0. That was why, in "old inflation," the bubbles must collide to fill the 
universe with energy; and the fact that this does not happen (because the 
bubbles grow only at the speed of light while the space between them grows 
superluminally) was fatal for that approach to infiation (Guth & Weinberg 
1983) .0 But now this defect is turned into a virtue by arranging to have 
a second burst of infiation inside the bubble, to drive the curvature back 
toward zero, i.e., ~* 1- By tuning the amount of this second period of 
infiation, it is possible to produce any desired value of (Sasaki, Tanaka, 
& Yamamoto 1995; Bucher, Goldhaber, & Turok 1995; Yamamoto, Sasaki, 
& Tanaka 1995). The old problem of too much inhomogeneity beyond 
the horizon producing too large a value of the quadrupole anisotropy is 
presumably solved because the interior of the bubble produced in the first 
infiation is very homogeneous. 

I personally regard this as an existence proof that infiationary models 
producing Qq ~ 0.3 (say) can be constructed which are not obviously wrong. 
But I do not regard such contrived models as being as theoretically attractive 
as the simpler models in which the universe after infiation is predicted 
to be fiat. (Somewhat simpler two-inflaton models giving fig < 1 have 
been constructed by Linde & Mezhlumian 1995.) Note also that if varying 
amounts of infiation are possible, much greater volume is occupied by the 
regions in which more inflation has occurred, i.e., where Oq ~ 1- But the 
significance of such arguments is uncertain, since no one knows whether 
volume is the appropriate measure to apply in calculating the probability of 
our horizon having any particular property. 

The spectra of density fluctuations produced in infiationary models with 
r^o < 1 tend to have a lot of power on very large scales. However, when 
such spectra are normalized to the COBE CMB anisotropy observations, 
the spherical harmonics with angular wavenumber £ « 8 have the most 
weight statistically, and all such models have similar normalization (Liddle 
et al. 1996a). 



t Although there have been attempts to revive Old Inflation within scenarios in which the 
inflation is slower so that the bubbles can collide, it remains to be seen whether any such 
Extended Inflation model can be sufliciently homogeneous to be entirely satisfactory. 



Dark Matter and Structure Format- on 



65 



1.6.7 Inflation Summary 



The key features of all inflation scenarios are a period of superluminal 
expansion, followed by ("re-")hcating which converts the energy stored in 
the inflaton field (for example) into the thermal energy of the hot big bang. 

Inflation is generic: it fits into many versions of particle physics, and 
it can even be made rather natural in modern supersymmetric theories as 
we have seen. The simplest models have inflated away all relics of any 
pre-inflationary era and result in a flat universe after inflation, i.e., Q = 1 
(or more generally + il\ = 1). Inflation also produces scalar (density) 
fluctuations that have a primordial spectrum 



where V is the inflaton potential and rip is the primordial spectral index, 
which is expected to be near unity (near-Zel'dovich spectrum). Inflation 
also produces tensor (gravity wave) fluctuations, with spectrum 



where the tensor spectral index ~ (1 — Up) in many models. 

The quantity (1 — rip) is often called the "tilt" of the spectrum; the larger 
the tilt, the more fluctuations on small spatial scales (corresponding to large 
k) are suppressed compared to those on larger scales. The scalar and tensor 
waves are generated by independent quantum fluctuations during inflation, 
and so their contributions to the CMB temperature fluctuations add in 
quadrature. The ratio of these contributions to the quadrupole anisotropy 
amplitude Q is often called T/S = Qj/Ql; thus the primordial scalar 
fluctuation power is decreased by the ratio + T/S) for the same COBE 
normalization, compared to the situation with no gravity waves (T = 0). In 
power-law inflation, T/S = 7(1 — rip). This is an approximate equality in 
other popular inflation models such as chaotic inflation with V^(^) = m?(p'^ 
or A0^. But note that the tensor wave amplitude is just the inflaton 
potential during inflation divided by the Planck mass, so the gravity wave 
contribution is negligible in theories like the supersymmetric model discussed 
above in which inflation occurs at an energy scale far below mpi. Because 
gravity waves just redshift after they come inside the horizon, the tensor 
contributions to CMB anisotropies corresponding to angular wavenumbers 
£ 2> 20, which came inside the horizon long ago, are strongly suppressed 
compared to those of scalar fluctuations. The indications from presently 
available data (Netterfield et al. 1996; cf. Tegmark 1996, and Silk's article 




(1.12) 




(1.13) 
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in this volume) are that the CMB amphtude is rather high for £ 200, 
approximately in agreement with the predictions of standard CDM with 
h ~ 0.5, rjfo ~ 0.1, and scalar spectral index Up = 1. This suggests that there 
is little room for gravity-wave contributions to the low-£ CMB anisotropics, 
i.e., that T/S <^ 1. Thus tests of inflation involving the gravity- wave 
spectrum will be very difficult. Fortunately, inflation can be tested with the 
data expected soon from the next generation of CMB experiments, since 
it makes very specific and discriminatory predictions regarding the relative 
locations of the acoustic peaks in the spectrum, for example the ratio of the 
first peak location to the spacing between the peaks ^i/A^ 0.7 — 0.9 (Hu 
& White 1996, Hu ct al. 1997). 

On the other hand, inflation is also Alice's restaurant where, according 
to the Arlo Guthrie song, "...you can get anything you want ... excepting 
Alice". It's not even clear what "Alice" you can't get from inflation. It was 
initially believed that inflation predicts a flat universe. But now we know 
that you 

• can get < 1 (with A = 0), as discussed in the previous section. 

• can make models consistent with super gravity and the sort of 
four-dimensional physics expected from superstrings, in which case one 
may expect that inflation occurs at a relatively low energy scale, which 
implies T/S" ~ 0, a low reheat temperature implying no production of 
topological defects and presumably requiring that baryosynthesis occur 
at the electroweak phase transition, and plenty of inflation implying that 

r^o !~ 1. 

• can alternatively get strings or other topological defects such as textures 
during or at the end of inflation (e.g. Hodges & Primack 1991) — which 
however probably requires tuning of the inflation and/or string model, 
for example to avoid a fractal pattern of structure-forming defects, which 
would conflict with the observed homogeneity of structure on very large 
scales. 

And in many versions of inflation, the most reasonable answer to the 
question "what happened before inflation" appears to be eternal inflation, 
which implies that in most of the meta-universe, exponentially far beyond 
our horizon, inflation never stopped. 



1.7 Comparing DM Models to Observations: ACDM vs. CHDM 



Dark Matter and Structure Format-on 



67 



1.7.1 Building a Cosmology: Overview 

An effort has been made to summarize the main issues in cosmological 



model-building in Fig. Here the choices of cosmological parameters, 

dark matter composition, and initial fluctuations that specify the model are 
shown at the top of the chart, and the types of data that each cosmological 
model must properly predict are shown in the boxes with shaded borders 
in the lower part of the chart. Of course, the chart only shows a few of the 
possibilities. Models in which structure arises from gravitational collapse 
of adiabatic inflationary fluctuations and in which most of the dark matter 
is cold are very predictive. Since such models have also been studied in 
greatest detail, this class of models will be the center of attention here. 

Perhaps the most decisive issue in model building is the value of the 
cosmological expansion rate, the Hubble parameter h. If /i ~ 0.7 as some 
observers still advocate, and the age of the universe Iq ^ 13 Gyr, then only 
low-f^o models can be consistent with general relativity.^ Depending on 
just how large h and to a-re, a positive cosmological constant may also be 
necessary for consistency with GR, since even in a universe with 17 ^ the 



age to Hq^ = 9.78/1-1 Qyj. (g^^ Yig. Thus, with A = and l^o ^ 



h < 0.75(13Gyr/to)- The upper limit on h is stronger, the larger VLq is: 
with A = and Jlo > 0.3, h < 0.61(13Gyr/to); with A = and l^o > 0.5, 
h < 0.57(13Gyr/to)- It has been argued above that the evidence strongly 
suggests that Q,q > 0.3, especially if the initial fluctuations were Gaussian; 
thus, if we assume values of h = 0.7 and to = 13 Gyr, we must include a 
positive cosmological constant. For definiteness, the specific choice shown 
is J^o + = 1) corresponding to the flat cosmology inspired by standard 
inflation. 

In such a ACDM model, one might initially try the Zel'dovich primordial 
fluctuation spectrum, i.e. Psik) = Agk'^p with Up = 1. However, this 
might not predict the observed abundance of clusters when the amplitude 
of the spectrum is adjusted to agree with the COBE data on large scales. 
If Oq > 0.3, then COBE-normalized ACDM predicts more rich clusters 
than are actually observed. In that case, it will be necessary to change 
the spectrum. The simplest way to do that is to add some "tilt" — i.e., 

t It is important to appreciate that the possible to — Hq (age-expansion rate) conflict goes to the 
heart of GR and does not depend ruLcosmological-model-dependent issues Uke the growth rate 



of fluctuations. As explained in § 1.2 besides GR itself the only other theoretical input needed 
is the cosmological principle: we do not live in the center of a spherical universe; any observer 
would see the same isotropy of the distant universe, as reflected in particular in the COBE 
observations. That is enough to imply the Friedmann-Robertson- Walker equations, which give 
the to — Ho connection. GR is not just a theory whose intrinsic beauty and great success in 
describing data on relatively small scales encourage us to extrapolate it to the scale of the 
entire observable universe. It is the only decent theory of gravity and cosmology that we have. 
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Fig. 1.3. Building a Cosmological Model. (This figure was inspired by similar 
flow-charts on inventing dark matter candidates, by David Weinberg and friends, 
and by Rocky Kolb.) 



assume that Up < 1. This adds one additional parameter. One can also 
consider more complicated "designer" primordial spectra with two or more 
parameters, which as I have mentioned can also be produced by inflation. 
In any case, it is then also necessary to check that the large-scale cluster and 

galaxy correlations on large scales are also in agreement with experiment. 
As we will see, this is indeed the case for typical ACDM models. In order 
to see whether such a model also correctly predicts the galaxy distribution 
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on smaller spatial scales ^ 10/i~ Mpc, on which the fluctuations in the 
number counts of galaxies Ng arc nonlinear — i.e., {SNg/Ng)rias ^1 — it is 
necessary to do N-body simulations. As we will discuss, these arc probably 
rather accurate in showing the distribution of dark matter on intermediate 
scales. But simulations are not entirely reliable on the scales of clusters, 
groups, and individual galaxies since even the best available simulations 
include only part of the complicated physics of galaxy formation, and omit 
or treat superficially crucial aspects STich as the feedback from supcrnovae. 
Thus one of the main limitations of simulations is "'galaxy identification" — 
locating the likely sites of galaxies in the simulations, and assigning them 
appropriate morphologies and luminosities. 

If /i ~ 0.5 and to ^ 13 Gyr, or if /i ~ 0.6 and to ^ 11 Gyr, then models 
with critical density, Q = 1, are allowed. Since the COBE- normalized 
CDM model greatly overproduces clusters, it will be necessary to make 
some modification to decrease the fluctuation power on cluster scales — for 
example, tilt the spectrum or change the assumed dark matter composition. 
As we have discussed, hot dark matter cannot preserve fluctuations on small 
scales, so adding a little hot dark matter to the mix of cold dark matter and 
baryons will indeed decrease the amount of cluster-scale power. A possible 
problem is that tilting or adding hot dark matter will also decrease the 
amount of power on small scales, which means that protogalaxies will form 
at lower redshift. So such models must be checked against data indicating 
the amount of small-scale structure at redshifts z > 3 — for example, against 
the abundance of neutral hydrogen in damped Lyman a absorption systems 
in quasar spectra, or the protogalaxies seen in emission at high redshift. 
Acceptable models must of course also fit the data on large and small-scale 
galaxy distributions. As we will see, = 1 COBE-normalized models with 
a mixture of Cold and Hot Dark Matter (CHDM) can do this if the hot 
fraction ~ 0.2. 

The ultimate test for all such cosmological models is whether they 
will agree with the CMB anisotropics on scales of a degree and below. 
Such data is just beginning to become available from ground-based and 
balloon-borne experiments, and continuing improvements in the techniques 
and instruments insure that the CMB data will become steadily more 
abundant and accurate. CMB maps of the whole sky must come from 
satellites, and it is great news for cosmology that NASA has approved 
the MAP satellite which is expected to be ready for launch by 2001, and 
that the European Space Agency is planning the even more ambitious 
COBRAS/SAMBA satellite, recently renamed Planck, to be launched a few 
years later. 
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Both sorts of models that have been discussed — 0, = 1 tilted CDM 
(tCDM) or CHDM, and Oq + ^^A = 1 ACDM— are simple, one-parameter 
modifications of the original standard CDM model. The astrophysics 
community has been encouraged by the great initial success of this theory 
in explaining the existence of galaxies and fitting galaxy and cluster data 
(BFPR, DEFW), and the fact that biased CDM only missed predicting 
the COBE observations by a factor of about 2. The other reason why the 
CDM-variant models have been studied in much more detail than other 
cosmological models is that they are so predictive: they predict the entire 
dark matter distribution in terms of only one or two model parameters 
(in addition to the usual cosmological parameters), unlike non-Gaussian 
models based on randomly located seeds, for example. Of course, despite the 
relatively good agreement between observations and the predictions of the 
best CDM variants, there is no guarantee that such models will ultimately 
be successful. 

Although the cosmic defect models (cosmic strings, textures) are in 
principle specified in terms of only a small number of parameters (in the case 
of cosmic strings, the string tension parameter Gfi plus perhaps a couple of 
parameters specifying aspects of the evolution of the string network), in 
practice it has not yet been possible for any group to work out the predicted 
galaxy distribution in such models. Most proponents of cosmic defect models 
have assumed an O = 1, Hq ~ 50 cosmology, but the chart refers instead to 
a cosmic defects option under Hq = 70. This is done because it would be 
worthwhile to work out a low-J7 case as well, since in defect models there is 
less motivation to assume the infiation-inspired fiat (f^o+^A = 1) cosmology. 

1.7.2 Lessons from Warm Dark Matter 

As has been said, the chart in Fig. only includes a few of the possibilities. 
But many possibilities that have been examined are not very promising. The 
problems with a pure Hot Dark Matter (HDM) adiabatic cosmology have 
already been mentioned. It will be instructive to look briefiy at Warm Dark 
Matter, to see that some variants of CDM have less success than others in 
fitting cosmological observations, and also because there is renewed interest 
in WDM. Although CHDM and WDM are similar in the sense that both 
are intermediate models between CDM and HDM, CHDM and WDM are 
quite different in their implications. The success of some but not other 
modifications of the original CDM scenario shows that more is required 
than merely adding another parameter. 

As explained above, WDM is a simple modification of HDM, obtained by 
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changing the assumed average number density n of the particles. In the usual 
HDM, the dark matter particles are neutrinos, each species of which has 
rii^ = 113 cm~'^, with a corresponding mass of m(z^) = Qi^po/rii, = Q,^92h'^ 
eV. In WDM, there is another parameter, m/mo, the ratio of the mass of 
the warm particle to the above neutrino mass; correspondingly, the number 
density of the warm particles is reduced by the inverse of this factor, so that 
their total contribution to the cosmological density is unchanged. It is true of 
both of the first WDM particle candidates, light gravitino and right-handed 
neutrino, that these particles interact much more weakly than neutrinos, 
decouple earlier from the hot big bang, and thus have diluted number density 
compared to neutrinos since they do not share in the entropy released by the 
subsequent annihilation of species such as quarks. This is analogous to the 
neutrinos themselves, which have lower number density today than photons 
because the neutrinos decouple before e"'"e~ annihilation (and also because 
they are fermions). 

In order to investigate the cosmological implications of any dark matter 
candidate, it is necessary to work out the gravitational clustering of 
these particles, first in linear theory, and then after the amplitude of the 
fluctuations grows into the nonlinear regime. Colombi, Dodelson, & Widrow 
(1996) did this for WDM, and Fig. |1.4| from their paper compares the square 
of the linear transfer functions T{k) for WDM and CHDM. The power 
spectrum P{k) of fluctuations is given by the quantity plotted times the 
assumed primordial power spectrum Pp{k), P{k) = Pp{k)T{k)'^ . The usual 
assumption regarding the primordial power spectrum is Pp{k) = Ak^'P, 
where the "tilt" equals 1 — Up, and the untilted, or Zel'dovich, spectrum 
corresponds to rip = 1. 

One often can study large scale structure just on the basis of such linear 
calculations, without the need to do computationally expensive simulations 
of the non-linear gravitational clustering. Such studies have shown that 
matching the observed cluster and galaxy correlations on scales of about 
20-30 h-^ Mpc in CDM-type theories requires that the "Excess Power" 
EP « 1.3, where 



and as usual cr{r) = {5p/ p){r) is the rms fluctuation amplitude in randomly 
placed spheres of radius r. The EP parameter was introduced by Wright et 
al. (1992), and Borgani et al. (1996) has shown that EP is related to the 
spectrum shape parameter T introduced by Efstathiou, Bond, and White 
(1992) (cf. Bardeen et al. 1986) by T 0.5(SP)-3-3. For CDM and the 



EP = 



o-(25 /i-iMpc)/o-(8 /j-^Mpc) 



(1.14) 



((j(25 /i-iMpc)/a(8 h-mvc))sCDM ' 
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Fig. 1.4. The square of the hnear transfer function T(fc) vs. wavenumber k = 
(27r)/A (in units of h Mpc"^) for (a) Warm Dark Matter (WDM), and (b) Mixed 
Dark Matter (MDM — CHDM with ^ 1 neutrino species). (From Colombi, 
Dodelson, & Widrow 1996, used by permission.) 



ACDM family of models, T = fi/i; for CHDM and other models, the formula 
just quoted is a useful generalization of the spectrum shape parameter since 
the cluster correlations do seem to be a function of this generalized F, as 
shown in Fig. 1.5. As this figure shows, F ~ 0.25 to match cluster correlation 
data. Peacock & Dodds (1994) have shown that F ~ 0.25 also is required to 
match large scale galaxy clustering data. This corresponds to EP ~ 1.25. 
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Fig. 1.5. The value of the J3 integral for sCDM and a number of ACDM and 
CHDM models evaluated at i? = 20/i~"'^Mpc is plotted against the value of the 
shape parameter T defined in the text. As usual, Js^R) — J^^cci'r)r^dr, where ^cc 
is the cluster correlation function. The horizontal dotted line is the J3 value for the 
Abell/ACO sample. The squares connected by the dashed line correspond from left 
to right to CHDM with = 1 neutrino species and ^l^ = 0.5, 0.3, 0.2, 0.1, and 
(sCDM) ; the square slightly below the dashed line corresponds to CHDM with Ni, = 
2 and fli, = 0.2; all these models have 17 = 1, = 0.5. and no tilt. The triangles 
correspond (1-to-r) to ACDM with (rJo,/i) = (0.3,0.7), (0.4,0.6) and (0.5,0.6). The 
two circles on the left correspond to CDM with h — 0.4 and (1-to-r) tilt (1 — Up) = 
0.1 and 0. These points and error bars are from a suite of truncated Zel'dovich 
approximation (TZA) simulations, checked by N-body simulations. (From Borgani 
et al. 1996.) 
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Since calculating a{r) is a simple matter of integrating the power spectrum 
times the top-hat window function, 

a^{r) = / P{k)W{kr)k'^dk (1.15) 



the linear calculations shown in Fig. 1.4 immediately allow determination of 



EP for WDM and CHDM. The results are shown in Fig. 1.6, in which the 
lower horizontal axis represents the values of the WDM parameter m/mo 
(with m/rriQ = 1 representing the HDM limit), and the upper horizontal 
axis represents the values of the CHDM parameter This figure shows 
that for WDM to give the required EP, the parameter value m/rriQ ~ 1.5 — 
2, while for CHDM the required value of the CHDM parameter is Ojy ~ 
0.3. But one can see from Fig. that in WDM with m/rriQ ^ 2, the 
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Fig. 1.6. Excess power EP in the two models discussed that interpolate between 
CDM and HDM. Solid curve shows EP as a function of the WDM parameter 
m/mo; note how quickly it becomes similar to CDM. Dashed curve shows how EP 
for Mixed Dark Matter (MDM — CHDM with N^^ = 1 neutrino species) depends on 
The observationally preferred value is EP « 0.25. (From Colombi, Dodelson, 
& Widrow 1996, used by permission.) 



spectrum lies a lot lower than the CDM spectrum at k ^ 0.3h~^ Mpc (length 
scales A 20 /i~^Mpc). This in turn implies that formation of galaxies, 
corresponding to the gravitational collapse of material in a region of size 
~ 1 Mpc, will be strongly suppressed compared to CDM. Thus WDM will 
not be able to accommodate simultaneously the distribution of clusters and 
galaxies. But CHDM will do much better — note how much lower T{k)'^ is at 
k > 0.3/1 Mpc"^ for WDM with m/mo = 2 than for CHDM with = 0.3. 
Actually, as we will discuss in more detail shortly, CHDM with = 0.3 
turns out, on more careful examination, to have several defects — too many 
intermediate-size voids, too few early protogalaxies. Lowering to about 
0.2, corresponding to a total neutrino mass of about 4.6(/i/0.5)~^ eV, in a 
model in which N^, = 2 neutrino species share this mass, fits all this data 
(PHKC95). 

Probably the only way to accommodate WDM in a viable cosmological 
model is as part of a mixture with hot dark matter (Malaney, Starkman, 
&: Widrow 1995), which might even arise naturally in a supersymmetric 
model (Borgani, Masiero, & Yamaguchi 1996) of the sort in which the 



Dark Matter and Structure Format-on 



75 



gravitino is the LSP (Dimopoulos et al. 1996). Cold plus "volatile" dark 
matter is a related possibility (Pierpaoli et al. 1996); in these models, the 
hot component arises from decay of a heavy unstable particle rather than 
decoupling of relativistic particles. 

There are many more parameters needed to describe the presently 
available data on the distribution of galaxies and clusters and their formation 
history than the few parameters needed to specify a CDM-type model. Thus 
it should not be surprising that at most a few CDM variant theories can fit 
all this data. Once it began to become clear that standard CDM was likely to 
have problems accounting for all the data, after the discovery of large-scale 
flows of galaxies was announced in early 1986 (Burstein et al. 1986), Jon 
Holtzman in his dissertation research worked out the linear theory for a 
wide variety of CDM variants (Holtzman 1989; cf. also Blumenthal, Dekel, 
& Primack 1988) so that we could see which ones would best fit the data 
(Primack & Holtzman 1992, Holtzman & Primack 1993; cf. Schaefer & 
Shafi 1993). The clear winners were CHDM with rjj, ~ 0.3 if /i ~ 0.5, and 
ACDM with ~ 0.2 if « 1. CHDM had first been advocated several 
years earlier (Bonometto & Valdarnini 1984, Dekel & Aarseth 1984, Fang 
et al. 1984, Shafi & Stecker 1984) but was not studied in detail until more 
recently (starting with Davis et al. 1992, Klypin et al. 1993). 



1.7.3 ACDM vs. CHDM — Linear Theory 

These two CDM variants were identified as the best bets in the COBE 
interpretation paper (Wright et al. 1992, largely based on Holtzman 
1989). In order to discuss them in more detail, it will be best to start by 



considering the rather complicated but very illuminating Fig. y_/^, showing 
COBE-normalized linear CHDM and ACDM power spectra P{k) compared 
with four observational estimates of P(A;).|] Panel (a) shows the 17 = 1 
CHDM models, and Panel (b) shows the ACDM models. The heavy solid 
curves in Panel (a) are for h = 0.5 and = 0.05. In the middle section of the 
figure, the highest of these curves represents the standard CDM model, and 
the lower ones standard CHDM (Ni, = 1) with ^l^, = 0.2 (higher) and 0.3; 
the medium-weight solid curves represent the corresponding CHDM models 
with two neutrinos equally sharing the same total neutrino mass {Ni, = 2). 
Note that the Ni, = 2 CHDM power spectra are significantly smaller than 
those for N,y = 1 for k ~ 0.04 — 0.4/i Mpc""*^; this arises because for N,y = 2 

f The normalization is actually according to the two-year COBE data, which is about 10% higher 
in amplitude than the final four-year COBE data (Gorski et al. 1996), but this relatively small 
diff'ercnce will not be important for our present purposes. 
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the neutrinos weigh half as much and correspondingly free stream over a 
longer distance. The result is that N,y = 2 COBE-normalized CHDM with 
~ 0.2 can simultaneously fit the abundance and correlations of clusters 
(PHKC95, cf. Borgani et al. 1996). The light solid curve is CDM with 
T = nh = 0.2. 

The "bow" superimposed on these curves represents the approximate 
"pivot point" (cf. Gorski et al. 1994) for COBE-normalized "tilted" models 
(i.e., with Hp 7^ 1). and the error bar there represents the la COBE 
normalization uncertainty. The window functions for various spherical 
harmonic coefficients ag, bulk velocities Vr, and ag are shown in the bottom 
part of this figure (see caption). The bow lies above the an window because 
the statistical weight of the COBE data is greatest for angular wavenumber 
£ ~ 11 (cosmic variance is greater for lower i, and the 7° resolution of the 
COBE DMR makes the uncertainty increase for higher £) . 

The upper section of Panel (a) reproduces the curves for sCDM (top), 
= 0.2 = 1 CHDM, and T = 0.2 (light) P{k), compared with 
several observational P{k) (see caption). Beware of comparing apples to 
oranges to bananas! Note that the only one of these observational data sets, 
that of Baugh Sz Efstathiou (1993, 1994) (squares) is the real-space P{k) 
reconstructed from the angular APM data; that of Peacock & Dodds (1994) 
(filled circles) is based on the redshift-space data with a bias-dependent 
and O-dependent correction for redshift distortions and a model-dependent 
(Peacock & Dodds 1996, Smith et al. 1997) correction for nonlinear 
evolution; the others are in redshift space. Also, the observations are of 
galaxies, which are likely to be a biased tracer of the dark matter, while 
the theoretical spectra are for the dark matter itself. Moreover, as will 
be discussed in more detail shortly, the real-space linear P{k) are only 
a good approximation to the true real-space P{k) for k ^ 0.2h Mpc"^; 
nonlinear gravitational clustering maJies the actual P{k) rise about an order 
of magnitude above the linear power spectrum for k ^ Ih Mpc^^. Thus 
one can see that COBE-normalized sCDM predicts a considerably higher 
P{k) than observations indicate. COBE-normalized F = 0.2 CDM predicts 
a power spectrum shape in better agreement with the data, but with a 
normalization that is too low. But the P{k) for = 0.2 CHDM, especially 
with = 2, is a pretty good fit both in shape and amplitude. The fact 
that the linear spectrum lies lower than the data for large k is good news 
for this model, since, as was just mentioned, nonlinear effects will increase 
the power there. 

The three heavy solid curves in Panel (b) represent the P{k) for ACDM 
with h = 0.8, fib = 0.02 for = 0.1 (top, for k = O.OOl/i Mpc"^), 0.2, and 
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Fig. 1.7. Fluctuation power spectra for COBE-DMR-normalizcd models: Panel (a) 
n = 1 CDM and CHDM models, Panel (b) f^o + f^A = 1 ACDM models. The 
theoretical spectra are discussed in the text. The data plotted for comparison is 
squares - real space P{k) from angular APM data (Baugh & Efstathiou 1993), 
filled circles - estimate of real space P{k) from redshift galaxy and cluster data 
(Peacock & Dodds 1994), pentagons IRAS 1.2 Jy redshift space P{k) (Fisher 
et al. 1993), open and filled triangles - CfA2 and SSRS2 redshift space P{k) (da 
Costa et al. 1994). At the bottom of each panel are plotted window functions 
for CMB anisotropy expansion coefficients ae (Panel (a): quadrupole a2, and an; 
Panel (b) left to right: 02 for Og = 0.1, 0.3, and 1), bulk flows Vr, and the rms 
mass fluctuation in a sphere of 8h~^ Mpc erg. (Prom Stompor, Gorski, & Banday 
1995, used by permission.) 
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0.3 (bottom). The lighter curves are for the same three values of plus 0.4 
(bottom) with h = 0.5, fi^ = 0.05 (the large wiggles in the latter reflect the 
effect of the acoustic oscillations with a relatively large fraction of baryons) . 
Dotted curves are for sCDM models with the same pair of h values. The 
observational P{k) are as in Panel (a). 

Note that the power increases at small k as JIq decreases, with opposite 
behavior at large k. Also, the COBE-normalized power spectra are 
unaffected by the value of h for small k, but increase with h for larger k 
(the fact that the light h = 0.2 curve in Panel (a) is lower than sCDM 
reflects the same trend). The fact that the data points lie lower than any of 
the ACDM models for k ^ 0.02 is worrisome for the success of ACDM, but 
it is too early to rule out these models on this basis since various effects such 
as sparse sampling can lead the current observational estimates of P{k) to 
be too low on large scales (Efstathiou 1996). A better measurement of P{k) 
on such large scales k ^ 10~'^h Mpc~^ will be one of the most important 
early outputs of the next-generation very large redshift surveys: the 2° 
field (2DF) survey at the Anglo- Australian Telescope, and the Sloan Digital 
Sky Survey (SDSS) using a dedicated 2.5 m telescope at the Apache Point 
Observatory in New Mexico. P{k) is much better determined for larger k 
by the presently available data, and the fact that the linear Qq = 0.2 and 
0.3 curves lie higher than many of the data points for larger k means that 
these h = 0.8 models will lie far above the data when nonlinear effects are 
taken into account. This means that, unless some physical process causes 
the galaxies to be much less clustered than the dark matter ("anti-biasing"), 
such models could be acceptable only with a considerable amount of tilt — 
but that can make the shape of the spectrum fit more poorly. 

1.7.4 Numerical Simulations to Probe Smaller Scales 

"Standard" n = 1 Cold Dark Matter (sCDM) with h « 0.5 and a 
near-Zel'dovich spectrum of primordial fluctuations (Blumenthal et al. 1984) 
until a few years ago seemed to many theorists to be the most attractive of 
all modern cosmological models. But although sCDM normalized to COBE 
nicely fits the amplitude of the large-scale flows of galaxies measured with 
galaxy peculiar velocity data (Dekel 1994), it does not flt the data on smaller 
scales: it predicts far too many clusters (White, Efstathiou, Sz Prenk 1993) 
and does not account for their large-scale correlations (e.g. Olivier et al. 
1993, Borgani et al. 1996), and the shape of the power spectrum P{k) is 
wrong (Baugh & Efstathiou 1994, Zaroubi et al. 1996). But as discussed 
above, variants of sCDM can do better. Here the focus is on CHDM and 
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ACDM. The linear matter power spectra for these two models are compared 



in Fig. 1.5 to the real-space galaxy power spectrum obtained from the 
two-dimensional APM galaxy power spectrum (Baugh & Efstathiou 1994), 
which in view of the uncertainties is not in serious disagreement with either 
model for 10^^ < k < Ih Mpc"^ The ACDM and CHDM models essentially 
bracket the range of power spectra in currently popular cosmological models 
that are variants of CDM. 
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Fig. 1.8. Power spectrum of dark matter for ACDM and CHDM models considered 
here, both normalized to COBE, compared to the APM galaxy real-space power 
spectrum. (ACDM: f^o = 0.3, = 0.7, h = 0.7, thus to = 13.4 Gy; CHDM: 
n ^ 1, ilu = 0.2 in Nu = 2 V species, h — 0.5, thus t^ = 13 Gy; both models fit 
cluster abundance with no tilt, i.e., Up — 1. (From Primack & Klypin 1996.) 



The Void Probability Function (VPF) is the probability Po{r) of finding 
no bright galaxy in a randomly placed sphere of radius r. It has been 
shown that CHDM with = 0.3 predicts a VPF larger than observations 
indicate (Ghigna et al. 1994), but newer results based on our = 0.2 
simulations in which the neutrino mass is shared equally between = 2 
neutrino species (PHKC95) show that the VPF for this model is in excellent 
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agreement with observations (Ghigna et al. 1996), as shown in Fig. I.E. 
However, our simulations (Klypin, Primack, & Holtzman 1996, hereafter 
KPH96) of COBE-normahzed ACDM with h = 0.7 and = 0.3 lead to a 
VPF that is too large to be compatible with a straightforward interpretation 
of the data. Acceptable ACDM models probably need to have f^o > 0.3 and 
/i < 0.7, as discussed further below. 
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Fig. 1.9. Void Probability Function Po(^) for (left panel) CHDM with h = 0.5 and 
= 0.2 in Ny — 2 species of neutrinos and (right panel) ACDM with h — 0.7 and 
J7o = 0.3. What is plotted here is difference between the actual VPF and that for 
a Poisson distribution, divided by V{R) = 47ri?^/3. Each plot shows also Po{R) for 
five typical different locations in the simulations (dotted lines) to give an indication 
of the sky variance. Data points are the VPF from the Perseus-Pisces Survey, with 
3(7 error bars; the VPF from the CfA2 survey is very similar. We have chosen 
the 5th for which the Pq of each model best approaches the observational data. In 
the top-left panel, the heavy "T" symbols at the bottom sets the boundary of the 
region where the signal is indistinguishable from Poissonian. They are obtained 
from the 3(7 scatters among measures for 50 different realizations of the Poissonian 
distribution in the same volume as our samples. (From Ghigna et al. 1996.) 



Another consequence of the reduced power in CHDM on small scales 
is that structure formation is more recent in CHDM than in ACDM. As 
discussed above (in § 1.4.7| ), this may conflict with observations of damped 



Lyman a systems in quasar spectra, and other observations of protogalaxies 
at high redshift, although the available evidence does not yet permit a 
clear decision on this (see below). While the original Vlj^ = 0.3 CHDM 
model (Davis, Summers, & Schlegel 1992, Klypin et al. 1993) certainly 
predicts far less neutral hydrogen in damped Lyman a systems (identified 
as protogalaxies with circular velocities Vc > 50kms~^) than is observed, 
as discussed already, lowering the hot fraction to fi^y ~ 0.2 dramatically 
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improves this (Klypin et al. 1995). Also, the evidence from preliminary data 
of a fall-off of the amount of neutral hydrogen in damped Lyman a systems 
for z ^ 3 (Storrie-Lombardi et al. 1996) is in accord with predictions of 
CHDM (Klypin et al. 1995). 

However, as for all 17 = 1 models, h ^ 0.55 implies to ^ 12 Gyr, which 
conflicts with the pre-Hipparcos age estimates from globular clusters. The 
only way to accommodate both large h and large to within the standard 
FRW framework of General Relativity, if in fact both h ^ 0.65 and Iq ^ 13 
Gyr, is to introduce a positive cosmological constant (A > 0). Low-Qq 
models with A = don't help much with to, and anyway are disfavored by 
the latest small-angle cosmic microwave anisotropy data (Netterfield et al. 
1997, Scott et al. 1996, Lineweaver & Barbosa 1997; cf. Ganga, Ratra, & 
Sugiyama 1996 for a contrary view). 

ACDM flat cosmological models with Qq = 1 — Q\ ~ 0.3, where = 
A/{3Hq), were discussed as an alternative to 17 = 1 CDM since the beginning 
of CDM (Blumenthal et al. 1984, Peebles 1984, Davis et al. 1985). They 
have been advocated more recently (e.g., Efstathiou, Sutherland, & Maddox 
1990; Kofman, Gnedin, & Bahcall 1993; Ostriker & Steinhardt 1995; Krauss 
& Turner 1995) both because they can solve the Ho—to problem and because 
they predict a larger fraction of baryons in galaxy clusters than 0, = 1 models 



(this is discussed in § 1.4.5 above) 



Early galaxy formation also is often considered to be a desirable feature 
of these models. But early galaxy formation implies that fluctuations on 
scales of a few Mpc spent more time in the nonlinear regime, as compared 
with CHDM models. As has been known for a long time, this results in 
excessive clustering on small scales. It has been found that a typical ACDM 
model with h = 0.7 and = 0.3, normalized to COBE on large scales 
(this fixes as ~ 1-1 for this model), is compatible with the number-density 
of galaxy clusters (Borgani et al. 1997), but predicts a power spectrum 
of galaxy clustering in real space that is much too high for wavenumbers 
k = (0.4 — l)/i/Mpc (KPH96). This conclusion holds if we assume either 
that galaxies trace the dark matter, or just that a region with higher density 
produces more galaxies than a region with lower density. One can see 
immediately from Fig. 1.7 and Fig. that there will be a problem with 
this ACDM model, since the APM power spectrum is approximately equal 
to the linear power spectrum at wavenumber k ~ 0.6h Mpc~^, so there is no 
room for the extra power that nonlinear evolution certainly produces on this 
scale — illustrated in Fig. |L1^ for ACDM and in Fig. [TTll for CHDM. The 
only way to reconcile the = 0.3 ACDM model considered here with the 
observed power spectrum is to assume that some mechanism causes strong 
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anti-biasing — i.e., that regions with high dark matter density produce fewer 
galaxies than regions with low density. While theoretically possible, this 
seems very unlikely; biasing rather than anti-biasing is expected, especially 
on small scales (e.g., Kauffmann, Nusser, &: Steinmetz 1997). Numerical 
hydro-l-N-body simulations that incorporate effects of UV radiation, star 
formation, and supernovae explosions (Yepes et al. 1997) do not show any 
antibias of luminous matter relative to the dark matter. 

Our motivation to investigate this particular ACDM model was to have 
Hq as large as might possibly be allowed in the ACDM class of models, which 
in turn forces 0,q to be rather small in order to have to ^ 13 Gyr. There is 
little room to lower the normalization of this ACDM model by tilting the 
primordial power spectrum Pp{k) = Ak^p (i.e., assuming Hp significantly 
smaller than the "Zel'dovich" value Hp = 1), since then the fit to data 
on intermediate scales will be unacceptable — e.g., the number density of 
clusters will be too small (KPH96). Tilted ACDM models with higher f^O) 
and therefore lower Hq for to }t 13 Gyr, appear to have a better hope of 
fitting the available data, based on comparing quasi-linear calculations to 
the data (KPH96, Liddle et al. 1996c). But all models with a cosmological 
constant A large enough to help significantly with the Hq — to problem are 
in trouble with the observations summarized above providing strong upper 
limits on A: gravitational lensing, HST number counts of elliptical galaxies, 
and especially the preliminary results from measurements using high-redshift 
Type la supernovae. 

It is instructive to compare the JIq = 0.3, h = 0.7 ACDM model that we 
have been considering with standard CDM and with CHDM. At k = 0.5h 
Mpc-\ Fi gures 5 and 6 of Klypin, Nolthenius, & Primack (1997) show that 
the r^jy = 0.3 CHDM spectrum and that of a biased CDM model with the 
same cg = 0.67 are both in good agreement with the values indicated for the 
power spectrum P{k) by the APM and CfA data, while the CDM spectrum 
with (Tg = 1 is higher by about a factor of two. As Fig. 1.11 shows, CHDM 
with = 0.2 in two neutrino species (PHKC95) also gives nonlinear P{k) 
consistent with the APM data. 



1.7.5 CHDM: Early Structure Troubles? 

Aside from the possibility mentioned at the outset that the Hubble constant 
is too large and the universe too old for any ^1 = 1 model to be viable, the 
main potential problem for CHDM appears to be forming enough structure 
at high redshift. Although, as mentioned above, the prediction of CHDM 
that the amount of gas in damped Lyman a systems is starting to decrease 
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Fig. 1.10. Comparison of the nonlinear power spectrum in the CIq = 0.3, h = 0.7 
ACDM model with observational results. Dots are results for the APM galaxy 
survey. Results for the real-space power spectrum for the CfA survey are shown 
as open circle (101/i~^ Mpc sample) and triangles (130/i~^ Mpc sample). Formal 
error bars for each of the surveys are smaller than the difference between the open 
and filled points, which should probably be regarded as a more realistic estimate 
of the range of uncertainty. The full curve represents the power spectrum of the 
dark matter. Lower limits on the power spectrum of galaxies predicted by the 
ACDM model are shown as the dashed curve (higher resolution ACDAl y simulation 
in KPH96) and the dot-dashed curve (lower resolution ACDM^. simulation). 

at high redshift z ^ 3 seems to be in accord with the available data, the 
large velocity spread of the associated metal-line systems may indicate that 
these systems are more massive than CHDM would predict (see e.g., (Lu 
et al. 1996, Wolfe 1996). Also, results from a recent CDM hydrodynamic 
simulation (Katz ct al. 1996) in which the amount of neutral hydrogen in 
protogalaxies seemed consistent with that observed in damped Lyman a 
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Fig. 1.11. Comparison of APM galaxy power spectrum (triangles) with nonlinear 
cold particle power spectrum from CHDM model considered in this paper (upper 
solid curve). The dotted curves are linear theory; upper curves are for z = 0, lower 
curves correspond to the higher redshift z = 9.9. (Prom Primack & Klypin 1996.) 



systems (DLAS) led the authors to speculate that CHDM models would 
produce less than enough DLAS (cf. Ma et al. 1997, Gardner et al. 
1997); however, since the regions identified as DLAS in these simulations 
were not actually resolved gravitationally, this will need to be addressed 
by higher resolution simulations for all the models considered before their 
arguments can be considered completely convincing. As mentioned above, 
fairly realistic simulations by Hachnclt, Stcinmctz, & Rauch (1997) arc able 
to account in rather impressive detail for the statistics characterizing the 
DLAS metal line systems measured by Prochaska &; Wolfe (1997), but 
Haehnelt et al. find that it is not large disks but rather clouds of gas in 
dark matter halos which account for most of these metal lines. 

Finally, Steidel et al. (1996) have found objects by their emitted 
light at redshifts z = 3 — 3.5 apparently with relatively high velocity 
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dispersions (indicated by the equivalent widths of absorption lines), which 
they tentatively identify as the progenitors of giant elliptical galaxies. 
Assuming that the indicated velocity dispersions arc indeed gravitational 
velocities, Mo &; Fukugita (1996, hereafter MF96) have argued that the 
abundance of these objects is higher than expected for the COBE-normalized 
n = 1 CDM-type models that can fit the low-redshift data, including 
CHDM, but in accord with predictions of the ACDM model considered 
here. (In more detail, the MF96 analysis disfavors CHDM with h = 0.5 
and rijy ^ 0.2 in a single species of neutrinos. They apparently would argue 
that this model is then in difficulty since it overproduces rich clusters — 
and if that problem were solved with a little tilt Up 0.9, the resulting 
decrease in fluctuation power on small scales would not lead to formation of 
enough early objects. However, if « 0.2 is shared between two species of 
neutrinos, the resulting model appears to be at least marginally consistent 
with both clusters and the Steidel objects even with the assumptions of 
MF96. The ACDM model with h = 0.7 consistent with the most restrictive 
MF96 assumptions has ^ 0.5, hence to ^ 12 Gyr. ACDM models having 
tilt and lower h, and therefore more consistent with the small-scale power 
constraint discussed above, may also be in trouble with the MF96 analysis.) 
But in addition to uncertainties about the actual velocity dispersion and 
physical size of the Steidel et al. objects, the conclusions of the MF96 
analysis can also be significantly weakened if the gravitational velocities 
of the observed baryons are systematically higher than the gravitational 
velocities in the surrounding dark matter halos, as is perhaps the case at 
low redshift for large spiral galaxies (Navarro, Frenk, &: White 1996), and 
even more so for elliptical galaxies which are largely self-gravitating stellar 
systems in their central regions. 

Given the irregular morphologies of the high-redshift objects seen in the 
Hubble Deep Field (van den Bergh et al. 1996) and other deep HST 
images, it seems more likely that they are mostly relatively low mass objects 
undergoing starbursts, possibly triggered by mergers, rather than galactic 
protospheroids (Lowenthal et al. 1996). Since the number density of 
the brightest of such objects may be more a function of the probability 
and duration of such starbursts rather than the nature of the underlying 
cosmological model, it may be more useful to use the star formation or 
metal injection rates (Madau et al. 1996) indicated by the total observed 
rest-frame ultraviolet light to constrain models (Somerville et al. 1997). The 
available data on the history of star formation (Gallego et al. 1996, Lilly ct 
al. 1996, Madau et al. 1996, Connolly et al. 1997) suggests that most of the 
stars and most of the metals observed formed relatively recently, after about 
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redshift z ^ 1; and that the total star formation rate at z ~ 3 is perhaps 
a factor of 3 lower than at z ~ 1, with yet another factor of ~ 3 falloff 
to z ~ 4 (although the rates at z ^ 3 could be higher if most of the star 
formation is in objects too faint to see). This is in accord with indications 
from damped Lyman a systems (Fall, Chariot, & Pei 1996) and expectations 
for n = 1 models such as CHDM, but perhaps not with the expectations for 
low-i^o models which have less growth of fluctuations at recent epochs, and 
therefore must form strT.ictT.irc earlier. BTit this niTTst be investigated TTsing 
more detailed modeling, including gas cooling and feedback from stars and 
supernovae (e.g., Kauffmann 1996, Somerville et al. 1997), before strong 
conclusions can be drawn. 

There is another sort of constraint from observed numbers of high-redshift 
protogalaxies that would appear to disfavor ACDM. The upper limit on 
the number of z ^ 4 objects in the Hubble Deep Field (which presumably 
correspond to smaller-mass galaxies than most of the Steidel objects) is far 
lower than the expectations in Iow-Qq models, especially with a positive 
cosmological constant, because of the large volume at high redshift in such 
cosmologies (Lanzetta et al. 1996). Thus evidence from high-redshift objects 
cuts both ways, and it is too early to tell whether high- or Iow-Oq models 
will ultimately be favored by such data. 

1.7.6 Advantages of Mixed CHDM Over Pure CD M Models 

There arc three basic reasons why a mixture of cold plus hot dark matter 
works better than pure CDM without any hot particles: (1) the power 
spectrum shape P{k) is a better fit to observations, (2) there are indications 
from observations for a more weakly clustering component of dark matter, 
and (3) a hot component may help avoid the too-dense central dark matter 
density in pure CDM dark matter halos. Each will be discussed in turn. 

(1) Spectrum shape. As explained in discussing WDM vs. CHDM 
above, the pure CDM spectrum P{k) does not fall fast enough on the 
large-A; side of its peak in order to fit indications from galaxy and cluster 
correlations and power spectra. The discussion there of "Excess Power" 
is a way of quantifying this. This is also related to the overproduction 
of clusters in pure CDM. The obvious way to prevent = 1 sCDM 
normalized to COBE from overproducing clusters is to tilt it a lot (the 
precise amount depending on how much of the COBE fluctuations are 
attributed to gravity waves, which can be increasingly important as the 
tilt is increased). But a constraint on CDM- type models that is likely to 
follow both from the high-z data just discussed and from the preliminary 
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indications on cosmic microwave anisotropics at and beyond the first acoustic 
peak from the Saskatoon experiment (Netterfield et al. 1997) is that viable 
models cannot have much tilt, since that would reduce too much both 
their small-scale power and the amount of small-angle CMB anisotropy. 
As already explained, by reducing the fluctuation power on cluster scales 
and below, COBE-normalized CHDM naturally fits both the CMB data 
and the cluster abundance without requiring much tilt. The need for tilt 
is further reduced if a high baryon fraction fti, ^0.1 is assumed (M. White 
et al. 1996), and this also boosts the predicted height of the first acoustic 
peak. No tilt is necessary for Cli, = 0.2 shared between Nn = 2 neutrino 
species with h = 0.5 and fib = 0.1. Increasing the Hubble parameter in 
COBE-normalized models increases the amount of small-scale power, so 
that if we raise the Hubble parameter to /i = 0.6 keeping = 0.2 and 
Ob = 0.1(0.5//i)^ = 0.069, then fitting the cluster abundance in this = 2 
model requires tilt 1 — Up ^ 0.1 with no gravity waves (i.e., T/S = 0; 
alternatively if T/S = 7(1 — Up) is assumed, about half as much tilt is 
needed, but the observational consequences are mostly very similar, with a 
little more small scale power). The fit to the small-angle CMB data is still 
good, and the predicted J7gas in damped Lyman a systems is a little higher 
than for the h = 0.5 case. The only obvious problem with h = 0.6 applies 
to any = 1 model — the universe is rather young: = 10.8 Gyr. But 
the revision of the globular cluster ages with the new Hipparcos data may 
permit this. 

(2) Need for a less-clustered component of dark matter. The fact 
that group and cluster mass estimates on scales of 1 /i^^Mpc typically 
give values for Q around 0.1-0.2 while larger-scale estimates give larger values 
around 0.3-1 (Dekel 1994) suggests that there is a component of dark matter 
that does not cluster on small scales as efficiently as cold dark matter is 
expected to do. In order to quantify this, the usual group M/L measurement 
of Oq on small scales has been performed in "observed" 0=1 simulations 
of both CDM and CHDM (Nolthenius, Klypin, & Primack 1997). We found 
that COBE-normalized O^ = 0.3 CHDM gives O^/l = 0.12-0.18 compared 
to ^M/L — 0-15 for the CfAl catalog analyzed exactly the same way, while 
for CDM ^M/L — 0.34 — 0.37, with the lower value corresponding to bias 
b = 1.5 and the higher value to 6 = 1 (still below the COBE normalization). 
Thus local measurements of the density in O = 1 simulations can give 
low values, but it helps to have a hot component to get values as low as 
observations indicate. We found that there are three reasons why this 
virial estimate of the mass in groups misses so much of the matter in 
the simulations: (1) only the mass within the mean harmonic radius r/^ 
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is measured by the virial estimate, but the dark matter halos of groups 
continue their roughly isothermal fallofF to at least 2r/j , increasing the total 
mass by about a factor of 3 in the CHDM simulations; (2) the velocities 
of the galaxies are biased by about 70% compared to the dark matter 
particles, which means that the true mass is higher by about another factor 
of 2; and (3) the groups typically lie along filaments and are significantly 
elongated, so the spherical virial estimator misses perhaps 30% of the mass 
for this reason. Visualizations of these simulations (Brodbcck et al. 1997) 
show clearly how extended the hot dark matter halos are. An analysis of 
clusters in CHDM found similar effects, and suggested that observations 
of the velocity distributions of galaxies around clusters might be able to 
discriminate between pure cold and mixed cold + hot models (Kofman et 
al. 1996). This is an area where more work needs to be done — but it 
will not be easy since it will probably be necessary to include stellar and 
supernova feedback in identifying galaxies in simulations, and to account 
properly for foreground and background galaxies in observations. 

(3) Preventing too dense centers of dark matter halos. Flores and 
Primack (1994) pointed out that dark matter density profiles with p{r) oc 
near the origin from high-resolution dissipationless CDM simulations 
(Dubinski &; Carlberg 1991; Warren et al. 1992; Crone, Evrard, & Richstone 
1994) are in serious conflict with data on dwarf spiral galaxies (cf. Moore 
1994), and in possible conflict with data on larger spirals (Flores et al. 1993) 
and on clusters (cf. Miralda-Escude 1995, Flores Sz Primack 1996). Navarro, 
Prenk, &; White (1996; cf. Cole & Lacey 1996) agree that rotation curves of 
small spiral galaxies such as DD0154 and DDO170 are strongly inconsistent 
with their universal dark matter profile PNFw{f) l/[r{r + o)^]. Navarro, 
Eke, & Prenk (1996) proposed a possible explanation for the discrepancy 
regarding dwarf spiral galaxies involving slow accretion followed by explosive 
ejection of baryonic matter from their cores, but it is implausible that 
such a process could be consistent with the observed regularities in dwarf 
spirals (Burkert 1995); in any case it will not work for low-surface-brightness 
galaxies. Work is in progress with Stephane Courteau, Sandra Faber, 
Ricardo Flores, and others to see whether the pnfw universal profile is 
consistent with data from high- and low-surface-brightness galaxies with 
moderate to large circular velocities, and with Klypin, Kravtsov, and Bullock 
to see whether higher resolution simulations for a wider variety of models 
continue to give Pnfw- The failure of earlier simulations to form cores 
as observed in dwarf spiral galaxies either is a clue to a property of dark 
matter that is not understood, or is telling us that the simulations were 
inadequate. It is important to discover whether this is a serious problem. 
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and whether inclusion of hot dark matter or of dissipation in the baryonic 
component of galaxies can resolve it. It is clear that including hot dark 
matter will decrease the central density of dark matter halos, both because 
the lower fluctuation power on small scales in such models will prevent the 
early collapse that produces the highest dark matter densities, and also 
because the hot particles cannot reach high densities because of the phase 
space constraint (Trcmaine & Gunn 1979, Kofman et al. 1996). But this 
may not be necessary, since even our simulations without any hot dark 
matter appear to be consistent with the rotation curves observed in dwarf 
irregular and low surface brightness galaxies (Kravtsov et al. 1997). 

1.7.7 Best Bet CDM-Type Models 

As said at the outset, the fact that the original CDM model did so well 
at predicting both the CMB anisotropics discovered by COBE and the 
distribution of galaxies makes it likely that a large fraction of the dark 
matter is cold — i.e., that one of the variants of the sCDM model might 
turn out to be right. Of these, CHDM is the best bet if turns out to be 
near unity and the Hubble parameter is not too large, while ACDM is the 
best bet if the Hubble parameter is too large to permit the universe to be 
older than its stars with 17 = 1. 

Both theories do seem less "natural" than sCDM, in that they are both 
hybrid theories. But although sCDM won the beauty contest, it doesn't fit 
the data. 

CHDM is just sCDM with some light neutrinos. After all, we know that 
neutrinos exist, and there is experimental evidence — admittedly not yet 
entirely convincing — that at least some of these neutrinos have mass, 
possibly in the few-eV range necessary for CHDM. Isn't it an unnatural 
coincidence to have three different sorts of matter — cold, hot, and baryonic 
— with contributions to the cosmological density that are within an order 
of magnitude of each other? Not necessarily. All of these varieties of 
matter may have acquired their mass from (super?)symmetry breaking 
associated with the electrowcak phase transition, and when we Tinderstand 
the nature of the physics that determines the masses and charges that are 
just adjustable parameters in the Standard Model of particle physics, we 
may also understand why 0^ and Jl^ are so close. In any case, CHDM 
is certainly not uglier than ACDM. 

In the ACDM class of models, the problem of too much power on small 
scales that has been discussed here at some length for ilg = 0.3 and h = 0.7 
ACDM implies either that there must be some physical mechanism that 
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produces strong, scale-dependent anti-biasing of the galaxies with respect 
to the dark matter, or else that higher and lower h arc preferred, 
with a significant amount of tilt to get the cluster abundance right and 
avoid too much small-scale power (KPH96). Higher Qq ^ 0.5 also is more 
consistent with the evidence summarized above against large and in favor 
of larger flo, especially in models such as ACDM with Gaussian primordial 
fluctuations. 

Among CHDM models, having N^, = 2 species share the neutrino mass 
gives a better fit to COBE, clusters, and small-scale data than N^, = 1, 
and moreover it appears to be favored by the available experimental data 
(PHKC95). But it remains to be seen whether CHDM models can fit the 
data on structure formation at high redshifts, and whether any models of the 
CDM type can fit all the data — the data on the values of the cosmological 
parameters, the data on the distribution and structure of galaxies at low and 
high redshifts, and the increasingly precise CMB anisotropy data. Reliable 
data is becoming available so rapidly now, thanks to the wonderful new 
ground and space-based instruments, that the next few years will be decisive. 

The fact that NASA and the European Space Agency plan to launch the 
COBE follow-up satellites MAP and COBRAS/SAMBA in the early years 
of the next decade, with ground and balloon-based detectors promising to 
provide precise data on CMB anisotropics even earlier, means that we are 
bound to know much more soon about the two key questions of modern 
cosmology: the nature of the dark matter and of the initial fluctuations. 
Meanwhile, many astrophysicists, including my colleagues and I, will be 
trying to answer these questions using data on galaxy distribution, evolution, 
and structure, in addition to the CMB data. And there is a good chance 
that in the next few years important inputs will come from particle physics 
experiments on dark matter candidate particles or the theories that lead to 
them, such as supersymmetry. 
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